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SUMMARY 


This  report  collates  and  reviews  the  results  of  an  eleven-year 
project  on  the  icing  of  fishing  trawlers. 

Initially  the  factors  causing  the  icing  of  ships  at  sea  are  presented 
and  some  aspects  of  the  physics  of  the  icing  process  are  discussed. 

Some  simple  icing  tests  on  cylinders  of  various  diameters  are 
described.  These  were  made  in  an  icing  wind  tunnel  and  demonstrated  the 
effects  of  air  temperature  and  cylinder  size  on  the  resulting  ice  formations. 
Later,  a  number  of  possible  methods  for  reducing  the  icing  hazard  were 
tested  in  the  wind  tunnel  and  at  an  outdoor  test  rig.  Of  the  methods  tested 
the  most  effective  was  an  inflatable  rubber  de-icing  blanket. 

Throughout  the  project,  icing  report  forms  were  distributed  to 
vessels  operating  out  of  Canadian  east  coast  ports.  The  data  obtained  by  this 
means  has  helped  to  establish  the  types  of  weather  systems  responsible 
for  icing  conditions,  the  geographical  extent  of  the  occurrence  of  icing,  and 
statistics  on  the  severity  of  icing  encounters  in  the  study  area. 

Appendices  present  simplified  formulae  for  the  droplet  collection 
efficiency  of  cylinders  and  rectangular  bodies,  the  derivation  of  an  analytical 
expression  for  the  rate  of  ice  build-up  on  fishing  trawlers,  and  a  bibliography 
of  ship  icing., 


(frangais  au  verso) 


RESUME 


Dans  le  present  rapport,  nous  compilons  et  analysons  les  resultats 
d’un  projet  de  onze  ans  sur  le  givrage  des  chalutiers. 

Dans  un  premier  temps,  nous  traitons  des  facteurs  a  l’origine  du 
givrage  des  navires  en  mer  et  de  certains  aspects  de  la  physique  du  processus 
de  givrage. 


Nous  decrivons  quelques  essais  simples  portant  sur  le  givrage  de 
cylindres  de  differents  diametres.  Ces  essais  ont  ete  menes  en  soufflerie  et 
ont  permis  dt  mettre  en  evidence  les  effets  de  la  temperature  de  l’air  et  des 
dimensions  du  cylindre  sur  ia  formation  de  la  glace  resultante.  Nous  avons 
par  la  suite  effectue,  en  soufflerie  et  sur  un  greement  d’essai  place  a  l’exte- 
rieur,  des  essais  sur  un  certain  nombre  de  methodes  susceptibles  de  reduire 
le  risque  de  givrage.  Parmi  ces  methodes,  la  plus  efficace  a  etc  l’utilisation 
d’une  couverture  de  degivrage  gonflable  en  caoutchouc. 

Tout  au  long  du  projet,  nous  avons  distribue,  aux  chalutiers  cir- 
culant  au  large  des  ports  de  la  cote  et  du  pays,  des  formules  de  releve  sur  le 
givrage.  Les  donnees  ainsi  obtenues  nous  ont  aide  a  determiner  les  types  de 
systemes  meteorologiques  qui  creent  des  conditions  de  givrage,  a  evaluer 
l’etendue  geographique  ou  se  produit  le  givrage,  et  a  produire  des  statistiques 
sur  l’intensite  du  givrage  observe  dans  la  zone  etudiee. 

Les  annexes  contiennent  des  formules  simplifiees  sur  le  coefficient 
de  captation  des  gouttelettes  des  cylindres  et  des  structures  rectangulaires, 
l’eiaboration  d’une  expression  analytique  donnant  le  taux  d’accumulation 
de  glace  sur  les  chalutiers  et  une  bibliographie  des  ouvrages  portant  sur  le 
givrage  des  navires. 
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TRAWLER  ICING 

A  COMPILATION  OE  WORK  DONE  AT  N.R.C. 


1.0  INTRODUCTION 

As  a  result  of  the  loss  of  three  Canadian  fishing  vessels  fiom  icing  in  preceding  winters,  a 
research  investigation  into  ship  icing  was  proposed  in  1968  jointly  by  the  then  Department  of  Transport, 
Ship  Safety  Branch,  and  the  National  Research  Council. 

This  project  included  a  preliminary  investigation  of  the  icing  characteristics  of  cylindrical 
objects,  and  some  simple  evaluation  tests  of  a  number  of  proposed  passive  de-icing  methods.  However, 
the  main  thrust  of  this  project  has  been  the  collection  of  data  on  icing  encounters  reported  by  ships 
operating  out  of  Canadian  east  coast  ports,  with  the  objective  of  determining  the  extent  and  severity 
of  icing  conditions  in  waters  off  the  Canadian  east  coast. 

This  report  consolidates  various  reports  and  papers  written  during  the  course  of  this  project 
which  has  now  been  terminated.  A  bibliography  is  presented  in  Appendix  C. 


2.0  FACTORS  CAUSING  THE  ICING  OF  SHIPS  AT  SEA 

A  perennial  hazard  faced  by  the  crews  of  vessels  operating  in  cold  northern  waters  during  the 
winter  months  is  that  of  icing.  This  applies  particularly  to  the  Canadian  east  coast  fisherman  for,  in 
winter,  conditions  conducive  to  icing  may  exist  at  any  time  after  leaving  port,  while  steaming  to  or 
from  the  fishing  grounds  or  while  actually  fishing;  in  fact,  his  whole  trip  is  likely  to  lie  within  a  zone  of 
potential  icing. 

The  icing  of  vessels  can  result  from  a  variety  of  causes  as  indicated  by  a  number  of  authors 
(Refs.  1  to  6)  as  follows: 

(a)  Supercooled  fog  (commonly  referred  to  as  arctic  frost  smoke  or  black  frost). 

(b)  Freezing  rain  or  freezing  drizzle. 

(c)  Falling  snow,  in  particular  wet  snow. 

(d)  Freezing  sea  spray  (vessel  and/or  wind  generated). 

Of  these,  the  freezing  of  sea  spray  on  a  vessel  is  the  most  dangerous  cause  of  icing.  The  other 
three  causes,  all  of  which  result  in  fresh  water  ice,  are  generally  not  considered  serious.  However, 
several  writers,  particularly  Borisenkov  (Ref.  3)  and  Shekhtman  (Ref.  8),  indicate  that  the  most  serious 
type  of  icing  occurs  when  one  or  more  of  these  three  lesser  causes  occurs  in  combination  with  freezing 
sea  spray.  On  the  other  hand,  it  has  been  suggested  by  Tabata  (Ref.  2)  that  snow  and  fog  should  not  he 
emphasized  in  connection  with  the  icing  of  ships. 

There  is  general  agreement  regarding  the  specific  meteorological  factors  that  result  in  the 
freezing  of  sea  spray  on  a  ship  (Refs.  3,  4,  7,  8).  The  air  temperature  must  he  helow  the  freezing 
temperature  of  the  sea  water  (about  -  1.8°C,  but  dependent  on  salinity),  the  sea  temperature  lower 
than  about  +6°C  (this  limit  is  less  well  defined  and  depends  on  other  factors,  particularly  air  tempera¬ 
ture),  and  the  wind  speed  greater  than  about  20  knots  (Beaufort  Force  5),  although  at  lower  wind 
speeds  some  light  icing  is  likely  to  occur. 

In  addition  to  these  factors,  which  apply  also  to  ocean  structures,  the  degree  of  icing  severity 
experienced  by  a  vessel  is  also  dependent  on  the  speed  and  heading  of  the  vessel  relative  to  the  wind  as 
well  as  on  the  size  and  lines  of  the  vessel  (including  freeboard,  amount  of  rigging  and  other  top  hamper, 
etc.). 
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Combined  with  the  relative  wind  speed  and  direction  factors,  the  wave  height  and  frequency 
also  have  a  bearing  on  the  severity  of  icing,  since  these  together  determine  the  amount  and  distribution 
of  water  and  spray  breaking  over  the  vessel.  Steep  waves  of  short  period  result  in  more  spray  than  long- 
period  waves  or  swell. 

It  is  apparent  from  all  this  that,  with  so  many  variables,  many  of  which  are  not  easily  meas¬ 
ured  or  characterized,  the  prediction  of  the  rate  of  icing  on  a  vessel  is  no  simple  task. 

3.0  THE  ICING  PROCESS 

The  formation  of  ice  on  a  body  results  when  water  impinges  on  the  body  in  a  sub-freezing 
environment  and  the  rate  of  loss  of  heat  from  the  body  is  such  that  some  or  all  of  the  water  is  frozen 
before  it  can  run  off  the  body.  Clearly,  the  colder  the  water,  the  greater  the  chance  that  this  will  occur. 
If  the  impinging  water  is  supercooled  (i.e.  its  temperature  is  below  its  equilibrium  freezing  point)  then 
it  is  a  certainty  that  ice  will  be  formed  providing  there  is  no  heat  input  to  the  process,  such  as  provided 
by  a  heated  body  or  by  the  kinetic  heating  of  high  speed  impact  between  the  water  and  the  body.  These 
cases  need  not  concern  us  in  a  discussion  of  the  icing  process  (although  a  body  might  be  heated  as  an 
icing  prevention  method). 

The  sea  water  may  be  several  degrees  above  its  equilibrium  freezing  point,  but  due  to  wave 
and  wind  action  and  to  the  passage  of  the  vessel,  water  is  detached  and  some  of  it  strikes  the  vessel. 

The  decks  may  be  washed  by  quantities  of  "green"  water  which  has  cooled  little  and  may  or  may  not 
create  much  ice.  As  height  increases,  the  quantity  of  water  becomes  less  and  its  degree  of  subdivision 
becomes  greater,  with  the  result  that  much  of  this  water  will  be  cooled  substantially,  and  the  smaller 
the  drop  size  the  more  nearly  will  its  temperature  approach  that  of  the  ambient  air.  These  drops  are 
likely  to  remain  in  a  supercooled  liquid  state*  during  their  flight  time  from  parent  sea  to  striking  the 
vessel's  structure.  In  general  then,  the  greater  the  height  above  the  sea.  the  lower  the  water  concentra¬ 
tion  will  be.  but  a  greater  proportion  of  that  water  will  freeze  and  form  ice. 

3.1  Drop  Deposition 

As  the  size  of  the  water  drops  decreases,  aerodynamic  drag  forces  increase  relative  to  the 
inertia  forces  acting  on  the  drops  and  there  is  a  tendency  for  the  trajectories  of  the  drops  to  more 
nearly  conform  to  the  streamlines  of  the  air  as  it  is  deflected  around  an  object.  The  larger  the  object, 
the  greater  the  deflection  for  a  given  drop  size  and  velocity  and  the  lower  the  proportion  of  drops  t  hat 
impinge  on  the  object.  Thus,  the  actual  impingement  or  catch  rate  is  less  than  the  potential  catch  rate 
assuming  straight  line  motion;  the  ratio  of  the  actual  to  potential  catch  is  known  as  the  catch  efficiency 
or  collection  efficiency.  Figure  1  shows  the  streamlines  of  air  around  a  cylinder  and  the  manner  in 
which  the  droplet  trajectories  are  deflected  from  a  straight  path.  Only  those  droplets  lying  between 
the  tangent  trajectories  impinge  on  the  cylinder,  thus  the  Collection  Efficiency,  Eln  ,  is  defined  by  the 

S 

ratio  —  . 

D 


As  the  drops  cool  below  their  equilibrium  freezing  temperature  they  do  not  immediately  freeze  but  exist  in  a  inctastjhlc  super¬ 
cooled  liquid  state.  Pure  water  can  support  a  maximum  supercooling  of  about  4(I0C  for  small  droplets  of  a  few  microns  diameter 
and  about  33 °C  tor  volumes  of  a  lew  cubic  centimetres  <  Ref.  9).  However,  in  nature  foreign  particles  in  the  water  act  as  freezing 
nuclei  causing  freezing  to  occur  at  temperatures  much  warmer  than  those  for  homogeneous  nuclcution.  nevertheless  significant 
supercooling  still  occurs,  the  most  likely  freezing  temperature  being  in  inverse  proportion  to  the  logarithm  of  droplet  size  I  he 
effect  of  salts  in  solution  on  the  other  hand  is  to  lower  the  nudeation  temperature  (Ref.  HU.  Ihe  effect  of  these  two  opposing 
phenomena  on  the  supercooling  and  nuclcution  of  sea  water  drops  is  unknown  but  is  obviously  ot  significance  to  the  study  of  sea 
spray  icing.  It  has  been  suggested  that  below  about  ~18T  the  spray  is  frozen  and  strikes  the  ship  as  small.  dr>  crystals  that  do  not 
adhere  (Ref.  I  l),  however,  this  is  contradicted  by  Shekhtman  (Kef.  8)  who  cites  instances  of  icing  occurring  at  air  temperatures 
down  to  -25°(\ 
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The  equations  of  motion  of  particles  in  a  stream  of  air  appear  to  have  been  first  developed 
by  Albrecht  (Ref.  12)  and  later  applied  by  Taylor  (Ref.  13)  to  droplets  causing  icing  of  aircraft.  Values 
of  Ein  are  not  easily  calculated;  however,  Langmuir  ( Ref.  14)  presented  droplet  trajectory  results  in 

graphical  form  for  simple  objects  such  as  spheres,  cylinders,  and  flat  ribbons  in  terms  of  two  non- 
dimensional  parameters.  The  use  of  graphical  results  of  this  sort  is  not  always  convenient,  particularly 
in  numerical  computations,  so  that  in  Appendix  A  (taken  from  Ref.  15)  empirical  expressions  are 
presented  that  permit  the  direct  computation  of  the  collection  efficiency  of  cylinders  and  flat  surfaces 
over  a  range  of  sizes  and  conditions  applicable  to  the  icing  of  ships  and  stationary  structures,  and  with 
sufficient  accuracy  in  view  of  present  uncertainties  with  regard  to  defining  the  spray  parameters. 


3.2  Heat  Balance 

When  the  supercooled  drops  strike  the  cold  obstacle,  it  acts  as  a  nucleating  agent  to  initiate 
freezing,  so  that,  depending  on  the  initial  supercooling  of  the  impinging  water  and  the  various  heat 
losses  and  gains  at  the  icing  surface  (Refs.  16,17)  some  or  all  of  the  impinging  water  freezes  in  its 
immediate  area  of  impingement.  The  proportion  that  does  actually  freeze  on  impact  is  called  the 
freezing  fraction,  n,  while  the  remaining  fraction  (1-n)  flows  away  from  its  impingement  zone  under 
the  influence  of  gravity  or  of  aerodynamic  forces;  some  of  this  run-off  will  flow  over  areas  of  greater 
heat  loss  where  further  freezing  will  occur  while  some  may  actually  flow  or  be  blown  completely  off 
the  surface  subject  to  impingement. 

The  freezing  fraction  and  hence  the  rale  of  icing  on  the  obstacle  is  given  by  the  solution  ot 
the  heat  balance  equation,  which,  if  only  the  primary  heat  transfer  processes  acting  at  the  icing  surface 
are  taken  into  consideration  (i.e.  ignoring  heat  transfer  due  to  radiation  and  kinetic  and  viscous  heating 
which  are  negligible  in  this  application ),  may  be  expressed  bv: 

q,  +  q*  +  q,  +  q,  +  q.,  =  o  n» 

where  q,  =  latent  heat  of  freezing  a  certain  fraction,  n,  of  the  impinging  water. 

qw  =  heating  (or  cooling)  of  the  impinging  water  to  the  equilibrium  surface 
temperature, 

qc  =  heat  transfer  by  convection  to  the  surrounding  air, 
qt,  =  heat  transfer  by  evaporation  to  the  surrounding  air. 

qa  =  heat  conducted  to  or  from  the  icing  surface  through  the  underlying  structure. 

The  fraction  of  the  impinging  water  that  freezes  decreases  both  as  the  temperature  increases 
and  as  the  water  concentration  increases.  The  excess  water  will  run  off  or  be  blown  off  the  surface,  and 
the  mass  of  ice  forming  on  the  obstacle  in  time,  t,  is 

Mj  =  nMw  =  nEm  VAwt  (2) 

where  n  =  freezing  fraction 

M„.  =  mass  of  water  impinging  on  obstacle  in  time,  t 
Kn)  =  collection  efficiency 

V  =  velocity  of  air  relative  to  obstacle 

A  =  projected  fronted  area  of  obstacle 

w  =  mass  of  water  droplets  in  unit  volume  of  air  (liquid  water  content ). 


As  the  ice  forms  on  the  obstacle,  the  shape  presented  to  the  flow  of  droplets  is  modified, 
and  the  projected  frontal  area  may  increase  in  size.  Thus  in  Equation  (2),  the  quantities  n,  Eni ,  and  A 

may  all  change  with  time,  so  that  the  instantaneous  icing  rate,  R,  =  dMj/dt,  also  becomes  a  function 
of  time. 

If,  in  Equation  (2),  A  is  considered  a  constant,  equal  to  the  frontal  area  of  the  ice-free 
obstacle,  then  it  is  possible  for  Em  and  even  the  Icing  Efficiency,  N  =  nEln  ,  to  assume  values  greater 

than  unity  if  the  frontal  area  of  the  iced  obstacle  exceeds  the  ice-free  area  A. 

The  above  discussion  has  looked  at  the  gross  or  overall  icing  rate  on  an  object,  and  where,  in 
Equation  (2),  the  factors  n  and  Em  are  the  effective  or  total  freezing  fractions  and  collection  efficiency 

for  the  ohjeet  as  a  whole.  No  information  is  forthcoming  from  such  considerations  on  the  local  growth 
rate  of  ice  on  different  parts  of  the  object,  and  hence  on  the  overall  shape  of  the  ice  deposit,  except  in  a 
rather  empirical  way  (Ref.  18). 

Determination  of  local  growth  rates  requires  a  knowledge  of  the  factors  relating  to  the  local 
convective  heat  transfer  and  to  the  local  collection  efficiency.  Such  factors  are  complex  and  are  far 
from  understood,  except  perhaps  for  the  simple  case  of  a  cylinder  for  which  an  attempt  (Ref.  19)  has 
been  made  to  numerically  model  the  local  icing  rate  and  so  predict  ice  shapes  on  a  cylinder. 

4.0  ICING  TESTS  ON  CYLINDERS 

As  an  introduction  to  the  National  Research  Council’s  programme  of  ship  icing  research,  a 
series  of  tests  was  made  in  an  icing  wind  tunnel  (1.37  m  square  test  section)  to  study  the  ice  formation 
and  rate  of  growth  on  cylinders  of  various  diameters.  These  tests  are  described  fully  in  Reference  17, 
and  their  salient  details  are  presented  here. 

Five  cylinders,  3.8  cm  (l'A  in.),  7.6  cm  (3  in.),  15.2  cm  (6  in.),  30.5  cm  (12  in.),  and  15.7  cm 
(18  in.)  in  diameter  were  used,  and  one-hour  tests  were  run  on  each  cylinder  in  both  a  horizontal  and 
vertical  orientation.  A  wind  speed  of  43  knots  (Beaufort  Force  9)  was  used  throughout  and  each  test 
was  repeated  at  two  temperatures  (about  -  14°C  and  about  -  7°C).  The  water  sprays  were  set  to  give 
a  water  concentration  in  the  air  of  3.2  g/m3  with  an  effective  drop  diameter  of  about  0.2  mm.  In  the 
light  of  the  paucity  of  quantitative  data  on  sea  spray  conditions,  these  values  were  chosen  rather  in¬ 
tuitively  within  the  capabilities  of  the  test  facility,  but  were  considered  to  result  in  fairly  representative 
icing  conditions.  One  area  where  conditions  where  clearly  not  representative  was  that  fresh  water  was 
used  rather  than  sea  water. 

Typical  ice  deposits  sure  shown  in  Figures  2  and  3,  and  representative  shapes  are  sketched  in 
Figure  4.  The  effect  of  greater  blow-off  and  run-off  of  water  at  higher  temperatures  is  evident  in  a 
smaller  ice  accretion  than  at  lower  temperature,  and  the  effect  of  gravity  results  in  an  asymmetrical 
formation  and  icicles  on  the  horizontal  cylinders,  and  a  variation  of  thickness  with  height  on  the 
vertical  cylinders. 

In  analyzing  the  results,  the  concept  of  Icing  Efficiency  was  used  (see  Section  3.2  above),  but 
in  the  form  of  a  one-hour  icing  efficiency,  defined  as: 


where 


mass  of  ice  formed  in  one  hour. 


The  ice  deposits  were  weighed  at  the  end  of  each  1-hour  run  and  the  1-hour  icing  efficiency 
determined.  In  spite  of  the  difference  in  the  shape  of  the  ice  accretions  on  horizontal  and  vertical 
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cylinders  (as  a  result  of  gravity  affecting  the  run-off),  no  significant  difference  in  icing  efficiency  was 
apparent.  It  was  therefore  possible  to  combine  the  results  from  both  orientations  as  in  Figures  5  and  6, 
which  show  the  effect  of  temperature  and  of  cylinder  diameter  on  the  1-hour  icing  efficiency. 


No  great  significance  should  be  read  into  the  actual  rates  of  icing  of  the  test  cylinders,  since 
their  relationship  to  rates  of  icing  at  sea  is  unknown.  However,  the  icing  efficiencies  derived  from  these 
tests,  and  presented  in  Figures  5  and  6,  are  considered  to  give  a  useful  guide  to  the  fraction  of  water  that 
actually  freezes  on  rails,  masts,  lines,  etc.  over  an  extended  period  of  time.  As  noted  in  Reference  20, 
“there  is  virtually  no  knowledge  as  to  the  rate  at  which  spray  strikes  the  exposed  surfaces  of  a  ship  or  of 
the  percentage  that  freezes”.  These  remarks  still  apply  some  20  years  later;  however,  these  results 
provide  some  indication  of  the  percentage  of  spray  that  may  freeze. 


The  manner  in  which  the  icing  efficiency  increases  with  decreasing  cylinder  size,  even  to 
values  greatly  in  exeess  of  100  percent  for  small  diameter  cylinders  (Fig.  6),  demonstrates  very  force¬ 
fully  the  absolute  necessity  for  keeping  the  number  of  rails,  stays,  etc.  to  a  minimum,  and  that,  where 
it  is  not  possible  to  dispense  with  them  completely,  one  larger  diameter  strut  should  be  made  to  do  the 
job  of  several  smaller  ones.  It  was  found  that  a  12-fold  increase  in  the  basic  cylinder  diameter  from 
3.8  cm  to  45.7  cm  resulted  in  only  a  3-fold  increase  in  the  mass  of  ice  collected. 

It  was  noticeable  that  the  ice  thickness  was  not  significantly  affected  by  the  cylinder  diameter, 
but  only  by  the  temperature.  On  the  other  hand,  the  width  of  the  ice  accretion  (dimension  A,  Fig.  4) 
was  affected  by  cylinder  diameter  as  well  as  temperature,  and  it  is  this  dependence  on  cylinder  diameter 
that  is  reflected  in  the  icing  efficiency.  Figure  7  demonstrates  the  effects  of  temperature  on  ice  thick¬ 
ness,  and  of  cylinder  diameter  on  the  ice  width  expressed  as  the  ratio  of  width  to  cylinder  diameter. 

The  fact  that  the  ice  thickness  appears  to  be  independent  of  cylinder  diameter  (a  relationship 
that  does  not  apply  when  all  the  water  freezes  in  its  area  of  impingement)  leads  to  the  speculation  that, 
for  sea  spray  ice,  the  rate  of  ice  growth  might  give  an  approximate  measure  of  the  spray  concentration 
providing  the  air  temperature  and  velocity  are  known.  Tests  at  other  water  concentrations  and  speeds 
are  required  to  ascertain  the  reliability  of  such  a  measure. 

The  main  findings  and  conclusions  from  these  tests  are: 

1.  The  dependence  of  the  icing  efficiency  on  the  cylinder  diameter  and  on  the  air  temperature 
was  demonstrated;  the  icing  efficiency  being  shown  to  increase  with  both  decreasing  dia¬ 
meter  and  decreasing  temperature. 

2.  It  was  shown,  for  the  conditions  of  the  tests,  that  the  ice  thickness  is  a  function  of  the 
temperature,  but  is  independent  of  the  cylinder  diameter. 

3.  The  width  of  the  ice  accretion  is  a  function  of  both  temperature  and  cylinder  diameter,  such 
that  for  small  diameter  cylinders  the  ice  width  may  grow  to  many  times  the  cylinder  diameter. 

4.  The  results  clearly  indicate  the  necessity  for  reducing  rigging,  spars,  etc.  to  a  minimum  as  one 
very  necessary  measure  for  diminishing  the  hazard  of  icing  at  sea. 


5.0  METHODS  OF  REDUCING  THE  ICING  HAZARD 

One  method  of  reducing  the  icing  hazard  was  demonstrated  in  the  tests  just  described,  that  is 
that  the  number  of  small  scale  collectors  such  as  rigging  and  railings,  and  top  hamper  generally,  should 
be  reduced  to  a  minimum  and  where  complete  elimination  is  not  feasible  one  larger  scale  component 
should  be  made  to  serve  the  function  of  several  small  scale  ones.  This  has  been  recognized  for  a  con¬ 
siderable  time  (Ref.  11),  but  that  it  has  not  always  been  fully  taken  to  heart  would  seem  to  be  evident 
from  some  of  the  discussion  in  the  final  report  of  a  British  inquiry  into  trawler  safety  (Ref.  21 ). 


Unless  a  vessel  can  completely  avoid  icing  conditions,  some  means  of  combatting  the  ice  un¬ 
necessary.  Methods  may  be  subdivided  into  two  categories,  de-icing  and  anti-icing.  Certain  de-icing 
methods  require  the  active  intervention  of  crew  members  on  deck:  these  include  such  methods  as  the 
use  of  steam  or  water  hoses  as  well  as  the  traditional  use  of  mallets  and  axes.  Conditions  on  deck  at 
these  times  are  particularly  hazardous,  and  become  completely  impossible  when  icing  conditions  an¬ 
tin'  most  severe  and  de-icing  is  the  most  necessary.  Also  such  methods  are  not  amenable  to  the  de-icing 
of  the  upper  parts  of  superstructure  and  rigging  where  the  added  weight  of  ice  contributes  most  to  tin- 
loss  of  stability.  Preferred  methods  therefore  should  preclude  the  use  of  crew  on  deck  and  so  must  be 
of  an  automatic  or  semi-automatic  nature.  Recent  proposals  or  developments  include  the  use  of  plastic 
foam  as  a  passive  de-icer  ( Refs.  22,  23),  inflatable  rubber  de-icer  blankets  ( Ref.  21 ),  and  an  electrical 
impulse  method  (Ref.  25).  The  cost  factor  seems  to  militate  against  the  extensive  use  of  electrical 
heating  (Ref.  26)  for  de-icing. 

Tests  at  NRC  (described  in  Ref.  27  and  summarized  below)  have  shown  the  remarkable 
effectiveness  of  the  inflatable  rubber  de-icer  on  both  cylindrical  and  flat  surfaces.  Little  is  known  about 
the  Russian  electro-impulse  system,  but  if  the  claims  made  for  it  are  valid,  it  might  well  lie  effective  for 
the  de-icing  of  large  flat  metal  surfaces. 

Anti-icing  methods,  in  which  ice  is  prevented  from  forming,  may  be  either  chemical  or 
thermal.  Chemical  methods  imply  basically  the  use  of  freezing  point  depressant  chemicals  applied  at 
the  surface  to  be  protected,  although  they  may  be  extended  to  include  materials  to  which  ice  has 
negligible  adhesive  strength.  The  major  problem  in  the  use  of  freezing  point  depressants  (i.e.  glycols 
and  alcohols)  is  that  of  uniformly  distributing  the  fluid  over  the  surface.  Materials  promoting  low  ice 
adhesion  are  theoretically  of  two  types,  fluid  and  solid,  the  fluid  type  being  typically  a  silicone  oil 
which  prevents  ice  from  adhering  to  the  underlying  solid  substrate;  such  a  method  although  resulting 
in  extremely  low  adhesion  (Refs.  28,  29,  30)  is  perhaps  only  practical  for  moderately  small  mechanical 
devices  on  which  an  oil  film  can  be  maintained  so  as  to  prevent  the  mechanism  from  being  jammed  by 
ice.  A  solid  material  to  which  ice  will  not  adhere  is  a  dream  not  yet  realized,  and  probably  unrealizable. 
Certain  fluoro-carbon  polymers  do  display  a  reduced  ice  adhesion  of  about  15  kN/m2  (20  psi)  in  shear, 
but  this  is  still  far  in  excess  of  the  approximate  1  kN/m2  necessary  to  promote  ready  ice  release. 

Certain  flexible  silicone  coatings  have  shown  low  values  (about  3.5  kN/m2 )  on  initial  ice  release,  but 
this  deteriorates  on  subsequent  releases  and,  in  addition,  problems  of  ultraviolet  degradation  and  low 
abrasion  resistance  currently  preclude  these  materials  from  practical  consideration. 

Thermal  anti-icing  methods  are  perhaps  rather  more  attractive  than  chemical  ones.  Large 
scale  electrical  heating  is  perhaps  prohibitively  expensive;  however,  for  small  items  it  may  provide 
certain  conveniences,  a  particular  case  in  point  being  the  use  of  an  electrically  conducting  transparent 
layer  deposited  on  glass  and  used  for  heating  bridge  windows  (Ref.  31,  Fig.  II).  It  seems  likely  that 
waste  heat  from  the  vessels'  main  engines  could  lie  employed  to  good  purpose,  either  to  heat  com¬ 
ponents  such  as  the  mast,  etc.  directly,  or  perhaps  indirectly  by  such  means  as  a  convective  thermo- 
syphon  as  described  by  Lock  (Ref.  32),  or  a  two-phase  thermosyphon  (Ref.  33),  the  action  of  which 
is  demonstrated  in  Figure  8.  The  two-phase  thermosyphon  has  the  advantage  of  higher  heat  transfer 
rates  and  of  remote  siting  of  the  heat  exchanger  and  boiler  section  as  long  as  the  vapour  transport  and 
liquid  return  are  unobstructed.  The  working  fluid  could  be  ammonia,  but  R21  would  be  more  suited 
to  the  temperatures  involved  since  the  pressures  in  the  system  would  then  be  dost'  to  atmospheric  and 
pressure  vessel  requirements  could  be  less  stringent.  Tests  at  NRC  have  demonstrated  the  anti-icing 
effectiveness  of  such  a  system  when  applied  to  simulated  masts  and  handrails. 

5.1  Tests  Made  at  NRC 

Evaluation  tests  (described  in  detail  in  Ref.  27)  were  made  on  a  number  of  methods  proposed 
for  reducing  the  hazard  of  superstructure  icing  on  ships  at  sea.  The  tests  were  made  under  simulated 
icing  conditions,  either  in  the  1.37  m  icing  wind  tunnel,  or  at  an  outdoor  test  site  using  natural  low 
temperatures  in  the  wintertime.  This  site  is  shown  in  Figure  9.  Ten  panels  1.2  m  X  0.9  m  were  attached 
to  a  vertical  frame  of  approximately  2.5  m  X  4.6m.  Various  methods  of  de-icing  were  mounted  on  or 
applied  to  the  panels.  Two  cables  (a  parallel  filament  rope  and  a  steel  cable)  were  strung  diagonally 


forward  and  down  from  the  two  upper  corners  of  the  frame.  To  augment  the  prevailing  wind  a  fan  was 
used  and  on  its  front  face  spray  nozzles  were  mounted  to  produce  the  icing  spray.  In  both  the  wind 
tunnel  and  outdoor  facility,  fresh  (tap)  water  was  used  to  produce  the  icing  conditions. 

The  methods  tested  were  all  de-icing  methods,  that  is  to  say  they  were  all  concerned  with  the 
removal  of  the  ice  after  it  had  formed  on  various  structures,  rather  than  with  the  prevention  of  its 
formation.  The  specific  methods  employed  were: 

1.  Pneumatic  de-icer  as  applied  to  a  mast 

2.  Plastic  foam  as  applied  to  a  mast 

3.  Plastic  foam  with  rubber  surface  as  applied  to  a  mast 

4.  Pneumatic  de-icer  as  applied  to  a  bulkhead,  such  as  a  bridge  front 

5.  Plastic  foam  as  applied  to  a  bulkhead 

6.  Plastic  foam  with  rubber  surface  as  applied  to  a  bulkhead 

7.  Paint  and  varnish  on  a  wooden  bulkhead 

8.  The  distribution  of  a  freezing  point  depressant  (ethylene  glycol)  over  a  bulkhead 

9.  Parallel  filament  rope  for  use  as  stays. 

5.1.1  Pneumatic  De-icer 

A  pneumatic  de-icer  is  essentially  a  tube  or  series  of  tubes  built  into  a  rubber  mat,  such  that 
when  the  tubes  are  inflated  with  air  they  expand  and  the  mat  surface  heaves  and  stretches,  and  in  so 
doing  breaks  down  the  adhesion  of  the  ice  to  the  mat  surface.  Such  pneumatic  de-icer  mats  have  long 
been  used  for  de-icing  the  leading  edges  of  aircraft  wings,  and  recently  their  application  has  been  ex¬ 
tended  to  marine  use.  In  this  application,  their  use  has  been  proposed  for  the  de-icing  of  masts  and 
stays,  bridge  fronts,  radar  scanners,  and  liferaft  stowages  (see,  for  instance,  Ref.  34  or  35). 

These  tests  constitute  an  independent  evaluation  of  their  effectiveness  in  de-icing  masts  and 
flat  surfaces.  The  test  specimen  was  a  short  length  of  mast  de-icer,  122  cm  long  and  91  cm  wide, 
suitable  for  attachment  to  a  30.5  cm  diameter  mast.  Details  are  given  in  Figure  10.  The  same  de-icer 
was  used  both  for  mast  de-icing  tests  in  the  wind  tunnel  and,  when  mounted  on  a  flat  surface,  for 
bridge  front  de-icing  tests  at  the  outdoor  test  stand. 

Six  tests  were  made  in  the  wind  tunnel  on  the  pneumatic  de-icer  installed  on  a  30.5  cm  dia¬ 
meter  cylinder,  137  cm  long,  to  represent  a  short  length  of  mast.  The  de-icer  was  not  bonded  to  the 
cylinder,  but  was  laced  onto  it  in  the  manner  shown  in  Figure  11.  The  cylinder  was  mounted  in  the 
tunnel  in  such  a  way  that  it  could  be  oscillated  about  its  normal  position  through  ±180°;  in  this  way 
ice  could  be  built  completely  around  the  mast,  a  situation  that  might  occur  should  a  vessel  change 
course  a  number  of  times  during  an  icing  encounter.  Thus  it  was  possible  to  see  whether  there  was  any 
change  in  the  de-icer’s  performance  between  complete  and  incomplete  encapsulation  of  the  mast  in  ice. 
Quite  low  temperatures  (-9  to  -23°C)  and  a  high  water  concentration  (3.2  g/nr' )  were  used  for  these 
tests  in  order  that  the  de-icer  be  subjected  to  a  severe  test. 

Only  two  of  the  wind  tunnel  mast  tests  will  be  described  here.  In  the  first,  at  23°(\  the 
mast  was  exposed  to  the  icing  condition  for  nearly  two  hours  during  which  time  it  was  rotated  *  180 
to  produce  an  ice  coating  uniformly  thick  around  the  whole  circumference.  The  radial  thickness  of 
the  ice  at  mid-height  was  about  13  cm  on  the  peaks  and  10  cm  in  the  valleys  of  the  rough  ice  formation 
(Fig.  12).  A  further  five  minutes  were  allowed  to  elapse  after  the  sprays  were  cut  off  to  allow  all  water 
on  the  surface  to  freeze  and  some  cooling  of  the  ice  to  occur.  One  inflation  of  the  de-icer  with  an  air 
pressure  of  100  kPa  (15  psi )  broke  this  thickness  of  fully  encapsulating  ice  cleanly  and  decisively,  as 
Figure  12  shows. 

The  second  test  to  be  described  was  at  a  high  temperature,  1 0°C,  and  for  a  duration  of 
only  45  minutes.  Because  of  the  higher  temperature,  the  water  flowed  on  the  ice  surface  to  a  greater 
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extent  than  on  the  earlier  test.  The  mast  was  not  rotated  for  this  test  and  so  a  rather  flat-faced  accretion 
resulted  (Fig.  13),  being  about  4  cm  thick  at  the  centre  front,  but  increasing  to  9  cm  thickness  at  45°  on 
either  side,  and  tapering  off  to  nothing  at  about  ±100°.  An  air  pressure  of  70  kPa  on  the  de-icer  shed  this 
ice  cleanly  with  a  single  inflation. 

Eight  tests  were  made  at  the  outdoor  test  site  on  the  pneumatic  de-icer  configured  as  a  bulk¬ 
head  de-icer.  The  de-icer  was  mounted  flat  on  one  of  the  1.2m  X  0.9m  panels  by  clamping  to  the  panel 
the  bars  that  passed  through  the  attachment  ferrules  down  each  side  of  the  de-icer.  The  de-icer  was  not 
bonded  to  the  panel  so  as  to  permit  the  required  deformation  when  the  inflation  air  pressure  was  applied. 

Because  of  the  variability  of  wind  direction,  intensity,  and  gustiness,  even  when  augmented  by 
the  fan,  it  was  not  possible  to  assess  an  effective  wind  velocity  for  each  run;  it  was,  however,  generally 
in  the  order  of  10  m/s  (20  kt).  For  this  reason,  and  because  throughout  these  tests  alterations  were  made 
in  the  water  spray  system  to  improve  the  spray  distribution  over  the  test  panels,  the  spray  concentration 
was  highly  variable  from  run  to  run.  Since  these  tests  were  very  qualitative  in  nature,  it  was  felt  that 
these  deficiencies  were  not  of  great  consequence  as  long  as  the  ice  formation  produced  appeared  reason¬ 
ably  characteristic  of  shipboard  icing. 

Again  only  representative  tests  will  be  described.  In  one,  the  de-icer  was  mounted  on  the 
upper  centre  panel.  About  10  cm  of  ice  was  allowed  to  form  at  a  temperature  of  -  15°C',  and  a  long 
overnight  cold  soak  at  -19°C  was  allowed  before  de-icing  was  attempted.  Thus  an  abnormally  severe 
test  was  given  the  de-icer,  requiring  more  than  one  inflation.  The  first  inflation  cracked  the  ice  on  the 
de-icer,  but  it  remained  anchored  around  the  edges  and  did  not  shed.  A  second  inflation  caused  a  small 
amount  of  ice  to  break  away  (Fig.  14b,  c)  but  not  until  the  third  inflation  was  the  bulk  of  the  ice  dis¬ 
lodged  (Fig.  14d).  Since  the  ice  was  anchored  at  the  edges,  and  particularly  at  the  top  where  the 
de-icer  tended  to  sag  away  from  the  steel  panel  behind,  there  seems  little  doubt,  had  the  de-icer  been 
much  larger  and  had  it  been  subjected  to  the  buffeting  experienced  at  sea,  that  the  ice  would  have 
shed  on  the  first  inflation. 

In  other  tests,  two  pneumatic  de-icers  were  used,  one  either  side  of  the  upper  centre  panel 
onto  which  was  bonded  rubber-faced  plastic  foam.  It  was  thought  that  the  simultaneous  inflation  of 
the  two  de-icers  might  also  remove  the  ice  from  the  panel  between.  Three  tests  were  made  with  this 
arrangement  at  de-icing  temperatures  of  -  12°C,  -8°C  and  -  3°C.  At  -  12°C,  ice  4  cm  thick  was  shed 
cleanly  from  one  de-icer  while  ice  on  the  other  de-icer  remained  anchored  along  the  top  and  sides, 
although  being  detached  from  the  active  surface  of  the  de-icer  (Fig.  15).  At  -8°C  ice  shed  cleanly  from 
both  de-icers,  but  only  at  -3°C  when  the  ice  was  rather  wet  and  rotten  was  the  desired  result  of  removing 
the  ice  from  the  intervening  rubber-covered  panel  achieved. 

5.1.2  Plastic  Foam 

As  a  result  of  Japanese  sea-going  tests,  mats  of  plastic  foam  applied  to  the  hull,  bridge  front, 
breakwaters,  deck,  etc.  are  reported  to  be  “effective  in  the  prevention  of  icing,  and  furthermore  make 
ice  removal  extremely  easy”  (Ref.  22). 

It  was  known,  before  Reference  22  was  available,  that  foamed  plastic  was  being  experimented 
with  in  Japan,  but  it  was  not  clear  exactly  how  it  functioned  as  an  ice  release  mechanism,  nor  the 
precise  materials  used.  However,  some  foamed  polyethylene  sheet  was  obtained  for  experimental  use. 
This  material  was  10  mm  thick  and  had  a  specific  gravity  of  0.05. 

This  plastic  foam  sheeting  was  tested  both  in  the  wind  tunnel  on  a  30  cm  diameter  steel 
cylinder  representing  a  mast,  and  on  flat  steel  plates  on  the  outdoor  test  stand.  In  both  instances  the 
foam  was  bonded  onto  the  steel  substrate  using  a  rubber-based  adhesive.  Tests  were  made  both  on  the 
bare  plastic  foam  and  with  a  0.8  mm  thick  sheet  of  neoprene  rubber  cemented  to  the  outer  surface. 

Two  test  runs  were  made  in  the  wind  tunnel  on  the  bare  plastic  foam  at  a  temperature  of 
about  14°C.  Other  conditions  were  as  for  the  pneumatic  de-icer.  The  cylinder  was  not  rotated  and 
the  ice  grew  to  a  thickness  of  3  cm  on  the  front  centreline  in  20  minutes  of  icing.  As  soon  as  the  tunnel 
was  shut  down  after  Run  1,  the  ice  was  assailed  with  a  wooden  club.  The  ice  was  rather  mushy,  and 
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although  it  tended  to  be  compressed  somewhat  by  clubbing,  it  was  removed  moderately  easily.  After 
Run  2,  20  minutes  were  allowed  for  the  ice  to  harden  before  removal  was  attempted.  This  time  the  ice 
broke  off  in  chunks  without  much  effort;  however,  no  large  area  of  ice  was  removed  in  one  piece, 
making  removal  rather  tedious. 

Next,  the  polyethylene  foam  was  covered  with  the  sheet  of  neoprene  rubber,  which  gave  a 
smoother  and  more  durable  surface  to  the  foam.  Two  more  icing  runs  were  made.  Thirty  minutes  were 
allowed  after  Run  3  for  the  ice  to  harden.  The  effort  required  to  remove  the  ice  was  judged  to  be  about 
the  same  as  for  Run  2.  Only  five  minutes  were  allowed  after  Run  4,  which  was  made  at  a  temperature 
of  -8°C,  and  the  ice  removal  was  judged  much  easier  than  after  Run  1. 

In  the  tests  at  the  outdoor  test  site  it  was  found  that  ice  could  only  be  removed  with  some 
difficulty  from  polyethylene  foam  bonded  to  a  flat  surface,  although  with  less  difficulty  than  from 
plain  painted  steel  panels.  The  compressible  foam  beneath  the  ice  layer  permitted  the  ice  to  be  cracked 
and  broken  with  a  club  or  hammer  more  readily  than  on  an  inflexible  steel  pane);  however,  the  broken 
pieces  of  ice  generally  remained  adhered  to  the  rather  rough,  accommodating  surface  of  the  foam,  and 
it  was  usually  necessary,  after  hammering,  to  pull  the  broken  pieces  off  the  foam  laboriously  by  hand. 

In  addition,  the  shattering  of  the  ice  frequently  cut  the  foam. 

Ice  removal  from  the  rubber-surfaced  foam  was  somewhat  less  difficult  and  showed  a  decided 
temperature  dependence.  After  a  run  at  about  -  19°C  it  was  noted  that  the  5  cm  thick  hard  porous  ice 
was  moderately  difficult  to  remove  using  a  wooden  club,  while  ice  on  the  steel  panels  was  practically 
impossible  to  break  and  remove  by  this  method.  Ice  was  finally  removed  by  a  sledge  hammer.  Kven 
then,  only  small  amounts  of  ice  could  be  broken  away  with  each  blow.  The  ice  was  judged  slightly 
easier  to  remove  at  about  16°C,  although  thin  ice  (about  1  cm)  just  splintered,  leaving  pieces  still 
adhered  to  the  ruhber. 

At  -8°C  ice  could  be  removed  without  too  much  trouble;  after  initially  breaking  it  with  the 
club,  the  ice  was  lifted  off  easily  by  hand.  At  -  3°C  the  ice  was  very  soft  and  wet  and  was  easily  removed 
from  all  panels. 


5. 1 .3  Paint  and  Varnish 

The  basic  finish  on  the  flat  steel  panels  at  the  outdoor  test  site  was  a  black  ant'-rust  paint. 
During  the  course  of  the  tests  one  of  the  panels  was  repainted  with  a  grey  deck  paint,  and  ;.  wooden 
panel  was  finished  half  in  the  grey  deck  paint  and  half  in  varnish. 

There  was  no  discernible  difference  in  the  removal  of  ice  from  any  of  the  steel  panels,  whether 
painted  with  the  black  anti-rust  paint  or  the  grey  deck  paint  they  were  all  equally  difficult.  Removal 
became  easier  as  the  temperature  increased,  and  little  difficulty  was  experienced  at  3JC  in  removing  ice 
from  any  of  the  panels. 

The  ice  was  not  quite  as  difficult  to  remove  from  the  wooden  panel,  but  this  was  thought  to 
be  due  to  its  lower  flexural  rigidity  compared  with  that  of  the  steel  panels.  If  was  judged  that  the  ice 
released  slightly  more  easily  from  the  grey  painted  half  of  the  panel  than  from  the  varnished  half, 
although  the  difference  was  so  slight  as  to  be  of  no  real  significance. 

5.1.4  Freezing  Point  Depressant 

Kthylene  glycol  was  distributed  over  one  area  of  the  vertical  steel  panel  wall  of  the  outdoor 
test  stand.  The  glycol  was  fed  under  gravity  from  a  tank  placed  80  cm  above  a  horizontal  distribution 
pipe  from  which  it  was  discharged  through  a  series  of  1  mm  diameter  holes  spaced  at  7.0  cm.  The 
discharge  holes  were  so  oriented  that  the  glycol  flowed  onto  the  panel  surface,  and  as  it  flowed  dow  n 
the  panel  surface  the  individual  streams  from  each  hole  tended  to  spread  and  join,  particularly  when 
the  surface  was  wet. 
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In  theory  the  freezing  point  depressant  should  mix  completely  with  the  impinging  water  and 
form  a  mixture  of  which  the  freezing  point  is  lower  than  the  ambient  temperature,  thus  resulting  in 
complete  freedom  from  icing.  In  practice  it  is  not  possible  to  ensure  ideal  distribution,  quantity,  or 
mixing  of  the  depressant.  Thus  ice  will  form,  but  this  ice  is  usually  infested  with  the  depressant  and  so 
is  mushy  and  has  poor  adhesion  to  the  surface  on  which  it  forms;  removal  is  therefore  easy. 

This  was  found  to  be  the  case.  In  one  run  at  ~8°C,  about  a  quarter  of  the  panel  was  free  of 
ice  at  the  end  of  the  test,  which  deposited  about  2.5  cm  of  ice  over  the  test  area.  Elsewhere  on  the  panel 
the  ice  was  just  standing  clear  of  the  panel  and  only  required  a  touch  to  collapse  it.  In  another  test  at 
10°C  glycol  was  fed  at  a  rate  of  1.1  litre/hr  to  six  panels,  i.e.  an  average  rate  of  1  litre/hr  per  6  m2 .  At 
the  end  of  the  run,  the  glycol-infested  ice  was  broken  away  from  the  panels  with  great  ease. 

5.1.5  Parallel  Filament  Rope  and  Steel  Cable 

A  parallel  filament  rope  was  compared  with  a  steel  cable.  This  plastic  rope  consists  of  a 
load-bearing  core  of  highly  compressed  parallel  filaments  of  polyethylene  terephthalate  (a  linear  poly¬ 
ester)  covered  by  a  black  polyethylene  sheath.  To  quote  the  manufacturer,  “A  twistless  rope  has  been 
evolved  which  is  light,  strong,  flexible,  maintenance-free,  and  most  important  —  highly  resistant  to 
stretch”.  Most  significant  from  the  icing  point  of  view  was  the  smooth  polyethylene  sheath  and  its 
flexibility,  which  suggested  the  possibility  of  low  adhesion  between  ice  and  rope. 

A  6  m  length  of  2-ton  rope  (11  mm  O.D.)  was  strung  from  the  top  right-hand  corner  of  the 
mounting  frame  for  the  wall  panels  to  the  front  left  corner  of  the  decking  (directions  as  viewed  from 
the  front),  while  for  comparison  afim  length  of  12.7  mm,  6X19  galvanized  steel  cable  was  strung 
from  the  top  left  corner  to  the  front  right  corner  (  Fig.  15). 

Because  of  the  difference  in  torsional  rigidity  between  the  parallel  filament  rope  anil  the 
steel  cable,  some  difference  in  the  manner  of  icing  was  observed.  The  rope,  because  of  its  lower  tor 
sional  stiffness,  tended  to  twist  under  the  weight  of  ice  forming  on  its  windward  side,  and  in  so  doing 
more  of  its  circumference  would  become  coated  with  ice  than  was  the  case  with  the  steel  cable  whose 
stiffness  precluded  significant  twisting  in  the  short  length  used.  This  twisting  effect  was  more  marked 
at  lower  temperatures  or  lower  water  concentrations  (Fig.  11)  under  which  conditions  more  of  the 
impinging  water  froze  on  the  windward  side  than  at  higher  temperatures  and  higher  water  concentra¬ 
tions,  when  most  of  the  ice  that  formed  on  the  cables  was  in  the  form  of  pendant  icicles  that  exerted 
little  torsional  moment  (Fig.  15). 

in  extreme  cases  complete  circumferential  encapsulation  of  the  plastic  rope  occurred  as  a 
result  of  the  twisting.  In  one  test  encapsulation  occurred  toward  the  lower  end  of  the  rope  so  that 
when  it  was  given  a  good  shake  by  hand  only  the  ice  on  the  upper  part  of  the  rope  was  dislodged. 

Since  the  rope  and  the  steel  cable  contacted  at  their  mid  points,  shaking  the  rope  also  shook  the  cable, 
with  the  result  that  most  of  the  ice  on  the  cable  was  also  shed.  The  ice  that  shed  had  left  the  rope 
completely  clean  and  free  of  ice  remnants,  indicating  true  adhesional  release,  while  small  ice  remnants 
still  adhered  to  the  steel  cable,  indicating  that  the  ice  had  failed  cohesively  rather  than  adhesively  at  the 
ice/metal  interface.  In  this  particular  test,  ice  on  the  lower  part  of  the  steel  cable  was  up  to  15  cm  thick, 
with  icicles  15  or  more  cm  long. 

It  was  found  on  this  and  other  occasions  that,  if  ice  at  the  lower  end  of  the  parallel  filament 
rope  was  removed  by  rapping  the  iced  rope  with  a  club,  encapsulating  ice  further  up  the  rope  would 
slide  down  and  could  in  turn  be  broken  off.  On  no  occasion  was  the  static  weight  of  the  ice  sufficient 
to  ove  -’ome  the  adhesive  force  between  the  ice  and  the  rope,  and  so  promote  a  self-protective  action. 


5.1.6  Discussion  and  Conclusions 

The  evaluation  of  these  various  methods  was  very  qualitative  and  the  test  conditions  may  not 
have  resembled  actual  sea  icing  conditions  with  any  degree  of  exactness.  However,  comparative  evalua¬ 
tions  are  thought  to  be  reasonably  valid.  The  use  of  fresh  water  instead  of  salt  water  is  thought  to  have 
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made  the  tests  somewhat  more  severe  by  virtue  of  the  greater  strength  of  fresh  water  ice.  Other 
mitigating  circumstances  at  sea  are  the  vibration  and  buffeting  experienced  in  rough  weather  —  condi¬ 
tions  that  could  not  be  simulated  on  the  simple  outdoor  test  facility. 


The  pneumatic  de-icer  proved  most  effective,  both  as  a  mast  de-icer  and  on  flat  surfaces. 

There  is  little  doubt  that  the  occasional  failure  to  shed  ice  on  the  first  inflation  of  the  de-icer  would  not 
have  occurred  under  the  buffeting  of  actual  sea-going  conditions  and  with  a  larger  active  area  than  used 
in  these  tests.  Disadvantages  of  this  otherwise  most  effective  method  of  de-icing  are:  1 )  cost,  2)  the 
likelihood  of  damage  to  the  de-icers  by  cutting  and  abrasion  if  used  in  or  adjacent  to  work  areas,  3 )  the 
danger  to  hands  working  on  deck  from  ice  falling  from  above,  and  4)  the  need  for  removal  of  the  shed 
ice  from  the  decks  if  reduction  of  the  total  ice  weight  is  to  be  achieved. 

On  the  other  hand,  pneumatic  de-icing  appears  to  be  an  ideal  method  of  de-icing  the  upper 
parts  of  masts,  etc.  to  reduce  the  large  heeling  moments  caused  by  ice  on  these  parts. 

The  use  of  ethylene  glycol  as  a  freezing  point  depressant  proved  effective  in  reducing  the 
strength  and  adhesion  of  ice  on  flat  surfaces.  It  is  probable  that,  under  the  buffeting  of  sea-going 
conditions,  a  fair  measure  of  anti-icing  performance  would  be  achieved  by  this  method  with  the  flow 
rates  used  in  these  tests.  This  method  of  icing  prevention  may  prove  useful  for  small  components  such 
as  radomes  and  inflatable  life  raft  containers.  Its  disadvantages  when  used  in  work  areas  are:  1 )  a  slip¬ 
pery  residue  could  be  left  on  the  decks  after  an  icing  encounter  if  not  completely  flushed  away  by  water, 
and  2)  a  danger  exists  of  contamination  of  the  fish  catch. 

The  various  coatings  and  finishes  tested  were  not  considered  very  successful.  The  most  suc¬ 
cessful  was  the  rubber -covered  plastic  foam.  This  was  because  the  deformable  foam  substrate  permitted 
breaking  of  the  ice  with  moderate  ease,  while  the  ice  exhibited  lower  adhesion  to  the  smooth  rubber 
surface  than  to  the  pebbled  surface  of  the  bare  polyethylene  foam.  The  rubber  also  provided  a  durable 
surface  to  the  otherwise  easily  damaged  foam.  No  case  occurred  of  ice  falling  off  the  foam  or  rubber- 
covered  foam  under  its  own  weight,  as  suggested  by  the  Japanese  work;  this  may  In'  due  partly  to  tin- 
small  panel  size,  partly  because  of  the  lack  of  vibration  and  buffeting,  and  also  perhaps  because  most 
of  the  tests  were  made  at  temperatures  well  below  the  freezing  point  with  fresh  water  instead  of  sea 
water. 


In  addition,  although  on  the  Japanese  vessels  of  steel  construction  the  mats  were  cemented 
on,  in  the  case  of  the  vessels  of  wooden  construction  the  mats  were  affixed  with  steel  straps  at  the  top, 
bottom,  and  sides  of  each  mat,  and  the  resulting  diaphragm  action  of  these  mats  under  the  action  of  the 
sea  undoubtedly  was  the  reason  that  “there  was  little  icing  on  the  mats,  and  ice  was  sprung  off  every 
time  there  was  any  kind  of  shock  or  shaking”. 

A  disadvantage  of  this  method  is  the  problem,  also  encountered  in  the  Japanese  tests,  of 
effectively  bonding  or  attaching  the  foam  and  rubber  blankets  to  flat  surfaces.  However,  little  difficulty 
is  seen  in  lacing  blankets  on  masts,  etc.,  much  in  the  manner  used  for  the  pneumatic  de-icer.  This  might 
prove  useful  in  work  areas  where  the  use  of  pneumatic  de-icers  may  be  prohibited  because  of  the  damage 
hazards. 


If  the  various  methods  are  rated  in  order  of  de-icing  effectiveness  or  ease  of  ice  removal,  the 
following  list  results: 

1.  Pneumatic  de-icer 

2.  Freezing  point  depressant 

3.  Rubber -surfaces  plastic  foam  on  steel  panel 

4.  Paint  or  varnish  on  wooden  panel 

5.  Bare  polyethylene  foam  on  steel  panel 

6.  Paint  on  steel  panel. 
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In  spite  of  the  obviously  lower  adhesive  strength  of  ice  to  the  polyethylene  sheath  of  the 
parallel  filament  rope  than  to  the  zinc  galvanizing  of  the  steel  cable,  little  difference  was  observed  in 
the  ease  of  removal  of  ice,  even  though  it  released  cleanly  from  the  rope,  while  cohesive  failure  of  the 
ice  resulted  in  residual  ice  fragments  remaining  on  the  steel  cable.  During  these  tests,  encapsulation 
occurred  on  the  plastic  rope  only,  because  of  its  lower  torsional  rigidity;  however,  in  stays  that  are 
more  vertical,  twisting  would  be  less  likely,  but  ice  encapsulation  might  be  more  likely  to  occur  as  a 
result  of  the  flow  of  water  around  the  cable.  Under  these  circumstances  difficulty  would  certainly  be 
encountered  in  removing  ice  from  steel  cables;  however,  the  tendency  for  ice  to  slide  down  the  rope 
once  ice  at  a  lower  level  had  been  removed  would  permit  the  removal  of  ice  from  otherwise  out-of-reach 
locations.  A  major  disadvantage  of  the  parallel  filament  rope  is  the  ease  with  which  it  can  be  cut,  thus 
limiting  the  locations  at  which  it  could  be  employed  aboard  ship. 


6.0  LOCAL  CONDITIONS  OFF  EASTERN  CANADA 

The  sea  area  of  interest  in  this  study  lies  approximately  between  the  latitudes  42°N  and  60°  N 
and  between  the  45°W  meridian  and  coast  of  North  America.  In  this  area  the  combination  of  strong 
winds  and  low  air  temperatures,  two  of  the  factors  required  for  icing,  are  provided  by  the  passage  of 
intense  cyclones.  Archibald  (Ref.  36)  states  that,  for  a  five-year  period  (1963-67)  for  January,  February 
and  March,  the  predominant  path  of  the  most  intense  storms  lies  right  along  the  Nova  Scotia  coast  and 
across  the  southern  part  of  Newfoundland,  with  a  somewhat  less  frequent  track  lying  parallel  and  south¬ 
eastward  of  this  predominant  track.  This  is  illustrated  in  Figure  16  which  is  reproduced  from  Archibald’s 
paper.  The  dots  with  concentric  circles  around  them  indicate  the  position  of  maximum  intensity  of  the 
centre  of  the  most  intense  storms,  i.e.  those  with  a  pressure  at  the  centre  of  between  956  and  972  mb 
and  with  a  very  steep  pressure  gradient.  Less  intense  storms,  with  pressures  of  between  972  and  980  mb, 
are  represented  by  solid  dots.  It  is  clear  that  many  of  these  storm  centres  will  produce  in  their  wake 
strong  northerly  or  northwesterly  winds  carrying  cold  Arctic  air  to  the  sea  area  of  interest.  The  result  is, 
therefore,  that  these  regions  are  particularly  susceptible  to  the  meteorological  conditions  most  conducive 
to  icing. 


The  oceanographic  conditions  of  interest  are  the  sea-surface  temperature,  the  sea  state  and,  to 
a  lesser  extent,  the  sea-surface  salinity. 

Monthly  mean  sea  surface  temperatures  are  presented  in  Reference  37.  Those  for  February, 
which  do  not  differ  very  significantly  from  those  for  January  and  March,  are  shown  in  Figure  17.  These 
show  cold  water,  with  temperatures  down  to  about  -1°C  adjacent  to  land  or  to  the  sea  ice  cover,  slowly 
increasing  in  temperature  with  distance  from  land  such  that  the  8°C  isotherm  lies  about  150  nautical 
miles  off  the  southern  tip  of  Nova  Scotia  (Cape  Sable),  some  250  nautical  miles  off  Cape  Race,  New¬ 
foundland,  and  considerably  more  removed  to  the  east  of  Newfoundland  and  Labrador. 

The  roughness  of  the  sea  depends  not  only  on  the  wind  strength,  but  also  on  the  fetch  and 
the  duration  of  the  wind  (Ref.  38).  This  is  illustrated  in  Table  I.  Some  sheltering  of  the  fishing  grounds 
off  Nova  Scotia  and  south  of  Newfoundland  from  northerly  winds  is  afforded  by  the  proximity  of  the 
land;  a  large  part  of  the  winter  fishing  is  done  within  100  nautical  miles  of  land,  resulting  in  significant 
reduction  in  potential  wave  height  for  winds  of  Force  6  and  above.  However,  a  large  part  of  the 
St.  Pierre  Bank  (45°  —  47°N,  55°  —  57°W)  and  the  Scatari  Bank  (46°N,  59°W)  is  not  protected  from 
the  north-westerly  winds  which  sweep  across  the  Gulf  of  St.  Lawrence  and  through  the  Cabot  Strait. 
Similarly  the  fetch  may  be  considered  unlimited  over  most  of  the  Grand  Banks  and  for  sea  areas  east  of 
Newfoundland. 

The  freezing  temperature  of  sea  water  is  a  function  of  its  salinity  and  may  he  represented  by 
the  expression: 


t  =-0.002  -  0.0524S  0.00006S2 


where  S 
and  t 


salinity  in  parts  per  thousand  (%<>), 
equilibrium  freezing  temperature  (°C). 
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This  expression  is  within  0.002°C  of  the  tabulated  values  of  the  freezing  (melting)  tempera¬ 
ture  of  sea  water  as  quoted  by  List  (Ref.  39)  for  salinities  from  0  to  40%, >. 

The  mean  isohalines  (Ref.  37)  for  the  three  winter  months  January  to  March  indicate  a  con¬ 
siderable  variation  in  the  area  of  interest,  from  about  30%<>  (parts  per  thousand)  in  the  Cabot  Strait  to 
about  35%n  at  the  southern  edge  of  the  Grand  Banks.  However,  this  variation  in  salinity  only  results  in 
a  variation  in  the  freezing  temperature  of  about  a  quarter  of  a  degree,  i.e.  from  about  -1.65°C  to  -1.9°C, 
as  Figure  18  demonstrates,  hardly  significant  enough  to  have  any  effect  on  the  onset  or  the  severity  of 
icing. 

These  considerations  lead  to  the  conclusion  that  conditions  conducive  to  frequent  and  severe 
icing  occur  in  the  area  under  consideration  during  the  winter  months. 

To  investigate  these  conditions  further,  to  compile  statistics  on  the  frequency  and  severity 
of  icing  encounters,  and  to  determine  actual  geographical  extents  of  icing  conditions,  data  have  been 
collected  by  means  of  voluntary  reports  using  a  standard  report  form.  The  next  section  will  discuss 
this  programme. 


7.0  ICING  REPORTING  PROGRAMME 

Icing  report  forms  were  distributed  jointly  by  the  Ship  Safety  Branch  of  the  Canadian  Coast 
Guard  and  the  Division  of  Mechanical  Engineering  of  the  National  Research  Council  of  Canada  to  all 
vessels,  but  primarily  fishing  vessels,  operating  out  of  Canadian  east  coast  ports.  Forms  were  distributed 
to  some  350  fishing  and  other  vessels  over  20m  in  length.  Of  these,  about  210  are  30m  or  longer,  and 
it  is  estimated  about  70%  may  be  at  sea  at  any  one  time.  Reporting  was  voluntary,  and  as  a  result  rather 
variable,  both  as  regards  compliance  and  the  quality  of  the  response;  however,  in  general,  both  the 
quantity  and  the  quality  of  the  responses  improved  with  time.  Reports  were  received  from  only  about 
15%  to  20%  of  the  vessels  over  30m  length,  while  on  any  one  day  no  more  than  about  10%  of  the  vessels 
actually  at  sea  made  out  reports,  even  when  icing  conditions  were  severe. 

The  report  forms  themselves  underwent  several  modifications  during  the  survey  period  in  an 
attempt  to  make  them  simpler  and  more  acceptable  to  the  user,  while  at  the  same  time  providing  as 
much  relevant  data  as  possible.  The  current  version  of  the  form  is  shown  in  Figure  19.  Earlier  versions 
were  presented  in  Refs.  15  and  40. 

Ships’  masters  were  urged  to  make  returns  even  when  no  icing  was  encountered  on  a  trip 
(i.e.  a  null  report)  so  that  a  better  appreciation  of  the  frequency  of  icing  encounters  could  be  gained. 

In  general  only  a  handful  of  masters  were  assiduous  in  this,  so  that  frequencies  of  occurrences  are  based 
only  on  a  proportion  of  the  vessels  responding,  the  number  of  which  varied  from  year  to  year,  as  did  the 
area  in  which  the  vessels  operated. 

Table  II  shows  the  number  of  vessels  responding  each  year  as  well  as  the  number  of  reports 
submitted. 

7. 1  Severity  Statistics 

Based  on  the  reports  from  all  vessels  for  seven  winters  and  on  the  reported  severity  of  icing 
(i.e.  the  subjective  assessment  given  by  the  ship’s  officer  making  the  report),  statistics  on  the  daily 
expectation  of  various  degrees  of  assessed  severity  for  each  of  the  winter  months  have  been  compiled. 
The  results  are  presented  in  Table  III.  For  each  of  the  months,  December  to  March,  and  for  the  com¬ 
plete  4-month  period,  1  December  to  31  March,  the  Table  gives  the  number  of  days  that  icing  of 
severity  light  or  greater,  moderate  or  greater,  or  heavy  was  reported.  The  yearly  average  number  of 
days  and  standard  deviation  from  this  average  is  shown.  The  great  variability  from  year  to  year  for  the 
months  of  December  and  March  is  evident  by  the  high  standard  deviations  (approximately  90%  of  the 
averages).  Much  greater  yearly  consistency  is  evident  for  the  months  of  January  and  February.  Based  on 
on  these  7  years  (1970-71,  1971-72,  1973-74,  1975-76,  1976-77,  1977-78,  and  1978-79)  the  daily 
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expectation  of  experiencing  icing  of  the  three  severity  categories  is  given;  for  example,  in  December 
more  than  trace  icing  may  be  expected  on  about  15 %  of  the  days  (i.e.  about  5  days)  while  in  February 
this  expectation  is  some  50^  (14  days).  Similarly,  for  the  4  month  season,  the  daily  expectation  is  36" 
(43  days)  for  light  icing  or  greater,  while  heavy  icing  may  be  expected  on  11  days  (9"). 

It  should  be  noted  that  the  data  of  Table  III  relate  to  icing  reported  from  throughout  the 
study  area.  Because  vessels  are  not  necessarily  in  every’  part  of  the  study  area  at  all  times,  these  statistics 
may  in  some  cases  rather  underestimate  the  frequency  of  occurrence. 

Frequencies  in  specific  areas  within  the  total  study  areas  may  be  expected  to  be  lower  than 
the  overall  frequencies  given  in  Table  III.  Two  specific  areas  have  been  selected  for  similar  analysis, 
namely  Divisions  3L  and  3P  of  the  area  designated  by  the  International  Commission  for  the  Northwest 
Atlantic  Fisheries  (ICNAF),  and  shown  in  Figure  22.  Table  IV  presents  the  results  lor  the  IC'XAF  Div¬ 
ision  3L  and  suggests  frequencies  of  about  one  third  of  those  for  the  total  study  area,  with  frequencies 
of  about  15"  (4  to  5  days  per  month)  for  icing  of  at  least  light  severity  during  the  months  of  January 
and  February.  In  February  the  frequency  drops  to  5.6"  for  moderate  or  greater  severity  with  large 
yearly  variability,  the  likely  cause  of  this  being  the  encroachment  of  sea  ice  into  this  area  to  a  greater 
extent  some  years  than  others.  The  results  for  the  ICNAF  Division  31*  (the  Sub-divisions  3 1'n  and  3l’s 
are  combined)  are  given  in  Table  V.  In  this  sea  area  the  frequencies  of  occurrence  are  seen  to  be  gener¬ 
ally  slightly  over  one  half  of  those  for  the  total  study  area,  except  for  the  month  of  December  when 
frequencies  of  less  than  2"  are  experienced.  The  higher  frequencies  for  December  seen  in  Table  III  for 
the  total  area  result  largely  from  icing  reports  from  more  northerly  sea  areas. 

The  foregoing  statistics  relate  to  the  subjective  icing  severity  assessments  as  made  by  the 
ships’  officers,  and  there  is  little  reason  to  question  these  assessments  in  view  of  the  great  experience 
these  men  have  of  icing  conditions  each  winter  and  of  the  ice  handling  qualities  of  their  respective 
vessels.  However,  these  statistics  relate  more  to  the  rate  of  icing  experienced  rather  than  to  the  total 
ice  accumulated  on  the  vessels.  Since  the  amount  of  accumulated  ice  is  the  significant  factor  as  far  as 
the  stability  and  handling  qualities  of  the  vessel  are  concerned,  and  is  of  special  significance  to  regulating 
bodies  concerned  with  the  safety  of  ships  at  sea,  statistics  based  on  ice  accumulation  have  boon  compiled. 

Since  most  nations  with  significant  deep  sea  fishing  fleets  subscribe  to  the  recommendations 
on  the  stability  of  fishing  vessels  of  the  Inter-Governmental  Maritime  Consultative  Organization  tlMCO). 
it  has  been  found  convenient  to  express  the  degree  of  icing  in  terms  of  the  'MCO  icing  allowance  crite¬ 
rion.* 


Based  on  seven  years  of  data  using  returns  from  only  those  vessels  for  which  a  relatively 
complete  history’  of  voyages  throughout  the  winter  season  was  available.  Table  VI  has  been  drawn  up 
and  updates  that  presented  in  Reference  12.  The  degree  of  icing  is  based  on  the  thicknesses  of  ice  on 
various  parts  of  the  vessel  as  reported  in  the  questionnaire  (Fig.  19).  and  converted  rather  arbitrarily  to 
multiples  or  sub-multiples  of  the  IMCO  criterion  using  the  points  system  illustrated  in  Figure  20.  The 
application  of  this  points  system  is  most  easily  demonstrated  by  the  use  of  an  example.  In  this  example 
it  will  be  assumed  that  the  following  ice  thicknesses  have  been  reported: 

a)  Diameter  of  ice  on  forward  rails  75-100  mm  (3-5  in.) 

b)  Diameter  of  ice  on  other  rails  75  100  mm  (3-5  in.) 


1  lie  Recommendation  on  Intact  Stability  ol  I  ishine  Vessels  t  Ret .  4  I  I  ot  the  IntcM  •overnmenta!  Maritime  (  onsulfativ  c  l  )iea ni/.i- 
tion  i  INK  Ol  lays  down  the  !  olio  wine  minimum  ice  allowance  in  calculating  a  vessel's  stabihtv 

a)  ke  in 2  <*n  exposed  weather  decks  ami  eaneways. 

h>  7  >  kf m;  h>r  the  projected  lateral  area  ot  ea*.h  suit  ot  the  vessel  above  the  water  plane 


(  n  ter  ton  a)  above  implies  an  aver  ace  ice  thickness  o(  about  b  min  oi  1 1  in  w  hile  hKoi  responds  to  an  av  erace  thk  knesx 
o|  roiiehlv  mm  or  J'X  m 
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c)  Diameter  of  ice  on  forward  stay 

d)  Thickness  on  main  deck 

e)  Thickness  on  boat  deck 

f)  Thickness  on  wheelhouse  font 

g)  Thickness  on  bulwarks 


>  100  mm  (>  5  in.) 
50-75  mm  (2-3  in.) 
25-50  mm  (1-2  in.) 
50-75  mm  (2-3  in.) 
25-50  mm  ( 1-2  in.) 


Thus,  using  Figure  20,  (a)  is  rated  1  points,  <b)  8  points,  (c)  8  points,  (d)  8  points,  (e)  1  points, 
(1)8  points  and  (g)  -1  )K)ints.  The  (Joints  total  is  therefore  -1-1,  which  falls  in  the  range  3(5-18,  thus  in¬ 
dicating  a  rating  equivalent  to  1'5  times  the  IMCO  icing  allowance. 

Fable  VI  shows  the  percentage  of  voyages  on  which  icing  of  the  indicated  degree  was  en¬ 
countered  by  the  assemblage  of  vessels.  If  icing  was  encountered  more  than  once  on  a  voyage,  only  the 
more  severe  encounter  was  used  in  these  staticties.  The  period  covered  is  a  typical  icing  season,  from 
roughly  the  middle  of  December  to  the  end  of  March  (approximately  100  days).  During  this  time  an 
individual  vessel  makes  some  7  or  8  trips  having  a  typical  duration  of  about  10  to  12  days. 

The  totals  for  the  seven  years  for  which  statistics  could  be  obtained  show  there  is  a  10% 
probability  of  icing  of  some  degree  being  encountered  on  a  voyage  to  the  fishing  grounds  during  the 
icing  season.  It  is  also  shown  that  on  8%  of  winter  voyages  the  degree  of  icing  encountered  will  exceed 
th(>  IMCO  icing  allowance,  while  2.5%  of  voyages  will  experience  icing  of  greater  than  twice  the  IMCO 
standard. 


7.2  Geographical  Distribution  of  Icing 

The  locations  of  icing  encounters  given  by  some  150  ship  icing  reports  are  shown  in  Figure1  21. 
They  indicate  that  icing  may  be  encountered  anywhere  within  the  continental  shelf  zone  of  Kastern 
Canada  as  well  as  the  adjacent  deeper  waters  of  the  Labrador  Sea.  The  very  large1  se>a  temperature  gra¬ 
dients  existing  off  the  Scotian  She'lf  and  the  Newfoundland  Grand  Banks  as  a  result  ejf  the  warm  Gulf 
Stream  provide  a  very  definite  limit  to  the1  southern  extern!  of  ic  ing  (Ref.  5). 

For  more  detailed  study,  the1  se>a  areas  in  which  icing  occurrences  have  been  reported  arc* 
conveniently  defined  in  terms  of  the  Divisions  designated  by  the1  International  Commission  fcjr  the1 
Northwest  Atlantic  Fisheries  (ICNAF).  Figure  22  is  a  map  showing  these  Divisions.  In  reporting  the 
are'as  of  icing  occurrence's  it  is  found  convenient  to  combine  Divisions  3l’n  and  31’s  into  a  single  Divi¬ 
sion  3I\  and  to  redefine  the  boundary  between  Divisions  IVn  and  IVs  by  moving  it  south  to  the  parallel 
1 110  N,  such  that  the  bulk  of  the  continental  shelf  within  -IV  is  in  the  mlefined  Division  -IVn. 

The  reports  suggest  that  there1  are  preferred  fishing  areas,  which  nevertheless  tend  to  vary 
somewhat  from  year  to  year  and  even  from  month  to  month,  partly  as  a  result  of  the  availability  of 
fish  stocks  and  also  as  a  re-sult  of  the  encroachment  of  se>a  ice,  both  southward  along  the  coasts  of 
Labrador  and  Newfoundland  and  within  the  Gulf  of  St.  I.awremee1.  Thus  it  is  not  possible  t o  interpret 
the  total  number  of  icing  reports  from  various  areas  as  indicating  the1  relative-  frceiucncy  and  or  severity 
of  icing  in  the1  respective  areas  at  any  givem  time.  I  lowe'ver,  by  taking  reports  of  icing  in  individual 
ICNAF  Divisions  fe>r  e*ach  e>f  the  winter  months  eeve'r  a  number  e>f  years  it  becomes  possible1  to  elete-r- 
mine-  whether  the-re1  is  a  likelihooel  ejf  cxpe-riencing  icing  in  a  givem  Division  and  memth.  Fight  ye-ars  erf 
elata  have1  been  analyzes!  by  Division  and  month  to  determine  whether  any  case-s  e>f  k  ing  of  severity 
greater  than  slight  (i.e1.  half  IMCO  erite-rion  eir  gre'ale'r)  have1  beam  re>port e‘e  1 .  The1  results  are1  prese-ntesl 
in  Table1  VII  anel  illustrated  as  pere-emtagecs  for  each  memth  in  Figure's  23  to  28. 

Since  a  require-el  e'emeiitiem  fe>r  se-rieius  icing  accoreling  t o  various  sejurces  (Re-fs.  3.  1,  7,  8)  is 
a  combination  e>f  air  te-mpe-rature  be'leiw  the1  fm-zing  point  ejf  se>a  wate’r  (about  1 ,8°C)  anel  se-a-surface- 
te'mpi’rature1  below  about  +5  (',  it  is  tei  lie1  expecte'el  that  such  a  combinatiem  e-ould  also  be1  use'll  to  ele- 
fine1  areas  of  e-x pe>ct e*<  1  icing.  Thus  in  Figure's  29  to  31  isotherms  of  10  F  (1. 1  C)  mean  monthly  se-a- 
surface1  ti'nipe-raturi1  as  given  in  Ri'fe'rime-e  37,  anil  isoph'ths  of  10%  pere-emtage-  fri'que-ne-y  of  air 
tempi-raturi-s  of  ()’C  or  le-ss  as  give-n  in  Ri'feremee1  13  are-  shown.  The  aetual  isotherms  of  both  sea  anel 
air  will  vary  throughout  e-ae'h  month  and  the-ir  monthly  means  will  vary  from  year  to  year:  however, 
comparing  Figure's  23  to  28  with  Figures  29  to  31  eloi-s  suggest,  that  thi'se1  monthly  means  takem  in 
combination  do  pre'ehet  the1  areas  whe-re1  icing  is  to  be-  expecteel  during  the1  month. 
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In  November  (Figs.  23  &  29)  no  icing  lias  been  reported  south  of  4H°N,  and  all  instances 
have  been  close  to  the  Labrador  or  Newfoundland  coasts,  observations  which  accord  with  the  sea  and 
air  temperature  criteria.  That  icing  has  not  been  reported  more  northerly  than  about  57  ' N  during 
November  is  doubtless  due  to  the  sparseness  of  Canadian  fishing  activity  in  these  waters:  nevertheless 
the  potential  for  serious  icing  in  the  Davis  Strait  area  exists.  Later  in  the  winter  the  southerly  advance 
of  sea  ice  tends  to  set  a  northern  limit  to  significant  icing  since  the  ice  cover  inhibits  wave  action  and 
hence  spray.  Light  to  moderate  icing,  however,  may  still  be  experienced  wherever  sufficient  open  water 
is  encountered.  Considerable  variability  in  ice  coverage  from  year  to  year  occurs. 

In  December  (Figs.  21  &  30)  the  possibility  of  encountering  icing  moves  further  south  with 
the  advance  of  the  cold  water  and  the  lowering  air  temperatures.  Probabilities  are  still  low  in  most 
areas,  and  no  encounters  have  been  reported  in  1CNAF  Divisions  3M,  3N  and  1  Vs.  The  single  report  of 
icing  in  Division  3-0  was  on  the  boundary  with  3L  and  3P.  so  that  the  13C  probability  indicated  may 
be  misleadingly  high. 

January  (Figs.  25  &  31)  shows  considerable  increase*  in  icing  probabilities,  with  the*  cold  water 
now  covering  all  of  the  Canadian  continental  shelf.  The  pattern  for  February  (Figs.  26  A  32 1  differs 
little  from  that  of  January  except  for  the  advance  of  sea  ice*  off  the  Labrador  Newfoundland  coast  and 
in  the  Ciulf  of  St.  Lawrence. 

March  (  Figs.  27  &  33)  sees  a  general  lessening  of  icing  probabilities  as  the  air  temperature 
begins  to  moderate,  while  only  isolated  cases  of  icing  have  been  reported  in  April  I  Figs.  2b  A  31 1. 

7.3  Meteorological  Conditions  Causing  Icing 

For  the  most  part  the  icing  reports  revealed  that  icing  occurrences  followed  closely  upon  the 
passage  of  well-developed  low  pressure  systems  through  the  area  (Ref.  51.  These  systems  bring  in  their 
wake  cold  Arctic  air  down  over  the  fishing  grounds  and  the  reports  indicate  that  they  influence  the*  area 
for  only  a  matter  of  a  day  or  two  at  a  time.  Figures  35  and  36  illustrate  such  systems.  Occasionally, 
however,  an  extensive  outbreak  of  cold  air  will  establish  itself  over  the  eastern  part  of  the  continent  for 
a  more  extended  period,  and  then  icing  conditions  may  persist  for  a  week  or  two  but  with  occasional 
lulls.  Systems  responsible  for  a  few  of  the  more  severe  icing  conditions  in  recent  years  are  illustrated  it' 
Figures  37  to  10. 

One  extended  period  of  icing  was  from  26  January  1972  to  10  February  1972.  and  Figure  37 
shows  the  weather  pattern  prevailing  at  1200  G.M.T.  on  27  January  1972.  Some  icing  was  reported  each 
day  throughout  this  period  with  severe  icing  being  experienced  on  26  to  28  January,  2  February.  6  Feb¬ 
ruary  and  9  February,  so  that  many  vessels  took  precautionary  action  such  as  dodging  into  wind  or 
sheltering  in  the  lee  of  land.  During  these  severe  periods  air  temperatures  as  low  as  1SJC  and  winds  of 
up  to  force  11  (60  knots)  were  reported,  while  the  sea  temperature  where  many  of  the  reporting  vessels 
were  operating  was  in  the  range  0  to  3  ('.  Such  conditions  can  result  in  very  severe  icing,  and  indeed 
this  is  borne  out  by  many  of  the  reports,  since  icing  rates  in  excess  of  15  cm  21  hr  were  evidenced  on 
several  occasions.  One  stern  trawler  of  630  gross  tons  reported  an  estimated  accumulation  of  200  tons 
of  ice  -  this  was  clearly  a  gross  overestimate,  but  ice  thicknesses  of  up  to  20  inches  (51  cm  )  occurred 
despite  several  de-icings  over  a  period  of  10  hours;  the  stability  of  the  vessel  was  reported  to  have  been 
noticeably  affected.  Another  stern  trawler  of  620  gross  tons  landed  at  Halifax  in  an  iced-up  condition: 
ice  from  several  locations  was  removed  and  weighed,  and  based  on  these  weights  and  the  thickness  of 
ice  at  various  positions  a  total  ice  weight  of  30.81  tons  was  arrived  at,  having  a  vertical  centre  of  gravity 
of  9.32  feet  (2.84  m)  above  deck  at  midships.  The  ice  allowance  for  this  vessel  in  its  stability  calcula¬ 
tions  is  30  tons  with  a  centroid  9.95  feet  (3.03  m)  above  the  level  of  the  deck  amidships.  The  greatest 
ice  thickness  on  this  vessel  was  12  inches  on  the  forward  house  front  (Fig.  11 ).  A  government  survey 
ship  of  65  m  overall  length  sailing  from  Newfoundland  to  Halifax  in  the  same  period  picked  up  an  esti¬ 
mated  150  tons  of  ice  and  reported  an  increased  roll  period;  Figure  42  shows  some  of  this  ice  build-up 
in  the  forecastle  area. 

These  occurrences  suggest  that,  compared  with  tne  IMCO  recommendations,  the  more  severe 
ice  allowance  imposed  by  Canada  (Ref.  44)  on  vessels  sailing  in  northern  or  eastern  Canadian  waters  is 
not  unrealistic  for  these  areas. 
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Another  extended  period  of  icing  was  the  12  days  from  9  February  1979  to  20  February 
1979,  with  paMicularly  severe  conditions  prevailing  on  14  and  15  February.  The  synoptic  situation  at 
0000  G.M.T.  on  14  February  is  shown  in  Figure  39  while  that  at  1800  G.M.T.  on  15  February  is  shown 
in  Figure  40.  The  more  southerly  low,  centred  at  about  40°N  55°W  at  0000  G.M.T.  on  the  14th,  deep¬ 
ened  and  moved  around  the  area  of  low  pressure  centred  just  off  the  Labrador  coast.  By  1800  G.M.T. 
on  the  15th,  the  latter  had  dissipated  off  the  northerly  tip  of  Newfoundland  while  the  former  had 
reached  its  greatest  intensity  and  was  centred  just  south  of  Greenland.  Throughout  this  period,  cold 
air  streamed  across  Labrador  and  the  sea  areas  off  Newfoundland  and  Nova  Scotia. 

7.4  Relation  Between  Reported  Severity  and  Temperature  and  Wind  Conditions 

To  demonstrate  the  effect  of  air  temperature  and  the  relative  wind  speed  on  the  reported 
icing  severity.  Figures  43  and  44  have  been  constructed  using  two  winters’  reports  (1971-72  and  1975- 
76),  for  vessels  steaming  at  5  knots  or  greater  (Fig.  43)  and  for  vessels  making  less  than  5  knots  ( Fig.  44  ). 
i.e.  when  riding  out  the  storm  (dodging)  or  when  fishing.  Although  best  fit  hyperbolae  have  been  drawn 
for  each  severity  level,  considerable  scatter  of  the  data  points  exists  and  no  clear-cut  transition  from  one 
severity  to  another  emerges;  nevertheless  the  trend  is  apparent.  The  scatter  may  be  attributed  to  the 
various  factors  enumerated  in  Section  2  above  as  well  as  to  the  subjective  nature  of  the  reported  icing 
severities. 


Figures  43  and  44  may  therefore  be  regarded  as  no  more  than  rough  rules  of  thumb  in  pre¬ 
dicting  icing  severity  if  only  the  air  temperature  and  wind  velocity  are  known.  Figure  43  suggests  that 
when  steaming,  heavy  icing  may  be  expected  if  the  product  of  the  relative  wind  speed  (in  knots)  and 
the  air  temperature  (in  °C  below  0°C)  is  greater  than  about  400,  while  when  dodging  or  fishing  (Fig.  44  ) 
the  criterion  for  heavy  icing  increases  to  a  product  in  excess  of  about  500.  Both  figures  suggest  that 
with  products  of  relative  wind  speed  and  air  temperature  less  than  200,  the  icing  severity  is  of  little  or 
no  significance. 

In  the  following  section,  an  attempt  is  made  to  quantitatively  predict  the  rate  of  ice  build-up 
on  a  vessel  using  an  analytical  approach. 


8.0  PREDICTION  OF  RATE  OF  ICING 

Certain  empirical  approaches  have  been  devised  to  determine  more  or  less  quantitatively  the 
severity  of  icing  from  the  prevailing  sea  and  atmospheric  conditions;  perhaps  the  best  known  of  these 
is  due  to  Mertins  (Ref.  45).  In  his  approach,  Mertins.  in  addition  to  air  temperature  and  wind  force, 
introduced  the  effect  of  sea-water  temperature  and  presented  his  results  in  the  form  of  four  charts 
( Fig.  45). 


More  recently,  Comiskey  (Ref.  46)  has  taken  Mertins’  four  charts  and  made  them  more  con¬ 
venient  to  use  by  combining  them  into  one  chart,  as  shown  in  Figure  46.* 

A  quantitative  rate  of  icing  may  be  approximated  from  both  the  Mertins  and  Comiskey  charts 
by  interpolation  using  the  range  of  icing  rates  given  for  each  icing  category. 

Applying  some  of  the  data  reported  by  means  of  our  icing  questionnaires,  the  Mertins/ 
Comiskey  charts  give  good  correlation  in  those  cases  when  the  vessels  are  steaming  (between  about 
5  knots  and  12  knots)  but  tend  to  underestimate  the  icing  rates  derived  from  the  reports  by  a  factor  of 
about  0.7.  In  cases  when  the  vessels  are  at  low  speed  (dodging  or  fishing)  correlation  is  not  so  good  and 
generally  the  charts  rather  seriously  overestimate  the  icing  rate.  Wind  speeds  relative  to  the  vessels  were 
used  rather  than  absolute  wind  speeds,  as  these  gave  somewhat,  improved  results. 
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Kachurin  ot  al  ( Ref.  47)  have  taken  an  analytical  approach  and  calculate  the  rate  of  icing  of 
a  cylindrical  rod  under  conditions  similar  to  the  spray  icing  of  small  fishing  boats,  and  relate  this  to  the 
ice  build-up  rate  observed  on  actual  ships.  This  approach  appears  to  be  useful,  but  as  formulated  by 
Kachurin  tends  to  be  rather  complex  inasmuch  as  he  introduces  time  dependent  effects. 

Such  sophistication  seems  to  be  rather  unjustified  in  view  of  the  many  uncertainties  in  meas- 
uting  or  estimating  the  various  atmospheric  and  ocean  parameters  as  well  as  a  representative  rate  of 
icing  for  the  vessel.  Therefore  in  Reference  42  an  expression  was  developed  for  estimating  the  rate  of 
icing  on  a  ship,  wherein  only  the  primary  heat  transfer  terms  were  considered,  and  where  in  addition 
the  ice  formation  process  was  considered  to  be  a  continuous,  quasi-steady  state  process,  notwithstanding 
that  the  water  arrives  at  the  icing  surfaces  in  a  rather  discontinuous  manner. 

In  Appendix  B  of  this  report  an  amended  version  of  the  analysis  given  in  Reference  42  is 
presented  wherein  the  drop  cooling  equation  is  modified  to  include  the  effect  of  evaporative  cooling, 
as  well  as  certain  other  minor  modifications. 

The  final  expression  for  the  rate  of  icing,  N(*,  derived  in  Appendix  B  is: 


N: 


*  =  0.00425  HwVkt  (ts  -  td)  +  0.0369  Vkt"  x  |  (ts  td)  +  17.3(e„  -  0.9  ea  )}  mm/h 


where 


Hw 

Vk, 


ea-es 


wave  height,  m 
relative  wind  speed,  kt 

equilibrium  temperature  of  the  icing  surface.  °C 
temperature  of  impinging  water  drops,  °C 
air  temperature,  °C 

saturation  vapour  pressure  of  moist  air  at  temperature  t  ,  t,  respectively,  kPa 


In  order  to  test  this  model  of  icing  under  conditions  of  sea  spray,  the  icing  rates  as  calculated 
by  this  equation  using  the  conditions  applicable  to  39  reported  icing  indieents  have  been  compared  with 
the  mean  icing  rates  derived  from  these  reports.  This  comparison  is  shown  in  Figure  47;  however,  before 
discussing  this  comparison,  some  discussion  of  the  data  on  which  it  is  based  is  necessary. 

8.1  Discussion  of  Observed  Data 

The  rfelevant  observed  and  derived  data  for  39  reported  icing  incidents  in  the  two  seasons 
1977-78  and  1978-79  are  presented  as  Table  VIII.  Only  data  from  fishing  trawlers  are  used;  those  from 
freighters,  tankers  and  government  vessels  of  various  types  have  not  been  included  because  of  the  dif¬ 
ferent  size  and  lines  of  these  vessels,  and  until  some  means  of  incorporating  these  two  variables  into 
the  analysis  is  developed,  each  class  of  vessel  must  be  treated  separately.  Even  among  the  fishing  trawl¬ 
ers  considerable  size  and  design  variation  exists  which  undoubtedly  results  in  part  in  some  of  the  scatter 
of  the  data. 

Scatter  may  also  be  caused  by  numerous  undertainties  in  the  data,  as  the  following  discussion 

indicates. 

The  reported  wind  speed  and  direction  are  combined  with  the  vessel  speed  and  heading  to 
derive  the  wind  speed  and  direction  relative  to  the  ship.  The  relative  wind  speed  is  entered  as  column  3 
of  Table  VIII.  Data  sets  (entries)  Nos.  1  to  32  are  for  relative  wind  directions  of  0  ±  45°,  while  entries 
33  to  39  are  for  those  cases  of  beam  or  quartering  winds  (i.e.  0  >  ±45°).  Some  uncertainty  exists  in 
interpreting  the  reported  data  since  in  many  cases  the  reported  wind  velocity  may  in  fact  be  the  relative 
velocity  rather  than  the  true  velocity,  and  it  is  not  always  made  clear  whether  the  ship  heading  is  mag¬ 
netic  or  true. 
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The  air  temperature  may  be  expected  to  be  fairly  reliable,  assuming  that  the  thermometer 
used  is  in  calibration  and  that  the  thermometer  bulb  is  not  exposed  to  the  water  spray.  Some  vessels 
are  evidently  not  equipped  with  thermometers,  and  so  data  from  them  could  not  be  used  except  where 
temperatures  could  be  inferred  from  other  vessels  operating  in  the  same  area  at  the  same  time. 

The  sea  water  temperature  entered  in  column  5  of  Table  VI 1 1  is  derived,  knowing  the  vessel 
location,  from  weekly  sea-surface  temperature  charts  compiled  by  the  Canadian  Forces  Meteorology/ 
Oceanography  Centre,  Halifax. 

The  wave  height  is  not  always  quantitatively  reported  and  in  most  instances  has  been  inferred 
from  a  qualitative  description  of  the  sea  state  using  the  following  scale  derived  from  the  State  of  Sea 
table  in  the  Trawlermen’s  Handbook  (Ref.  48). 


Description 

Inferred  Wave 
Height  —  m 

slight 

0.6  -  1.2 

moderate 

1.2  -  2.5 

rough 

2.5  -  4 

very  rough 

4  -  6 

high 

6  -  9 

very  high 

9  -  14 

If,  as  stated  in  Appendix  B,  the  liquid  water  content  of  the  spray  is  primarily  a  function  of 
wave  height,  then  this  is  an  important  parameter  in  the  analysis,  but  one  which  is  subject  to  considerable 
uncertainty  because  of  the  necessity  of  interpreting  the  qualitative  descriptions  given  in  the  reports. 
Relationships  relating  wave  height  to  wind  speed  exist,  but  these  usually  apply  to  winds  of  unlimited 
duration  and  fetch,  conditions  which  do  not  apply  to  the  majority  of  the  fishing  trawler  icing  events 
because  of  their  proximity  to  land  or  to  an  ice  cover. 

The  salinity  of  the  sea  water  in  the  location  of  an  icing  report  has  been  derived  from  the 
Mean  Sea  Surface  Salinity  charts  given  in  the  Oceanographic  Atlas  (Ref.  37).  Because  changes  in  salinity 
within  the  range  to  be  expected  in  this  area  (i.e.  30 %«  <  s  <  357(<>)  do  not  greatly  affect  the  icing  rate, 
secular  variations  from  the  salinity  values  given  in  the  charts  are  of  little  significance. 

Perhaps  the  greatest  uncertainty  arises  in  the  derivation  of  icing  rate  (column  8).  A  mean  or 
representative  ice  thickness  was  assessed  from  the  thicknesses  reported  fur  the  various  on-board  locations 
indicated  in  the  Icing  Report  (Fig.  19).  In  assessing  the  mean  thickness,  somewhat  less  weight  was  ac¬ 
corded  the  ice  thickness  on  the  forward  rails  and  stay,  since  these  thicknesses  were  not  necessarily  rep¬ 
resentative  of  the  vessel  as  a  whole.  The  reporting  of  duration  of  icing  occasionally  left  some  uncertainty, 
particularly  when  icing  was  intermittent,  in  which  case  what  should  have  been  reported  as  separate  events 
were  sometimes  reported  as  one  extended  event.  Sometimes,  too,  it  was  not  made  clear  whether  some 
de-icing  action  that  could  have  affected  the  reported  ice  thicknesses,  had  been  taken  during  the  icing 
encounter. 

Other  undertainties  arise  because  of  changing  conditions  during  the  encounter.  Often  such 
changes  (in  wind,  temperature,  etc.)  are  not  indicated.  Clearly,  the  longer  the  reported  encounter,  the 
more  uncertainties  there  are,  so  that  no  encounter  of  duration  greater  than  24  hours  has  been  used  in 
Table  VIII. 

8.2  Discussion  of  Icing  Rate  Results 

Employing  Equation  (B. 31c)  of  Appendix  Band  the  iterative  techniques  described  (Section 
B.1.7),  the  data  of  columns  2  to  7  of  Table  VIII  were  used  to  obtain  calculated  rates  of  icing  (column  9 
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of  Table  VIII).  These  are  compared  in  Figure  47  with  the  rates  of  icing  as  derived  from  the  icing 
reports  (column  8  of  Table  VIII).  Also  shown  in  the  figure  is  the  line  passing  through  the  origin  that 
gives  the  best  least  squares  fit  between  the  two  sets  of  values.  The  correlation  coefficient  is  r  =  0.55. 

That  there  is  so  much  scatter  in  Figure  47  is  perhaps  not  unexpected  in  view  of  the  numerous 
uncertainties  in  the  raw  data  discussed  in  sub-section  8.1.  It  is  encouraging  in  fact  that  the  calculated 
rate  is  generally  within  an  order  of  magnitude  of  the  assessed  rate,  which  gives  a  certain  degree  of  con¬ 
fidence  in  the  calculations;  however,  a  more  complete  test  of  the  theory  cannot  be  made  until  more 
consistent  and  reliable  data  is  available.  That  some  of  the  present  data  is  “bad”  may  be  seen  by  com¬ 
paring  some  of  the  entries  in  Table  VIII,  c.f.  entries  1  and  15,  25  and  33,  7  and  15,  10  and  21,  etc. 
where  lower  reported  icing  rates  correspond  to  clearly  more  severe  conditions. 

Figure  47  does  not  suggest  that  a  beam  wind  results  in  a  significantly  different  rate  of  icing 
than  a  bow  wind,  nor  can  any  consistent  effect  of  vessel  length  be  observed. 

The  rate  of  icing  Equation  (B.31c)  may  be  used  to  show  the  relative  effects  of  the  various 
parameters  contributing  to  the  heat  balance.  This  is  illustrated  in  Figure  48  where  one  parameter  of 
a  set  of  standard  conditions  has  been  varied  in  turn  throughout  its  likely  range.  Tin*  standard  conditions 
chosen  were: 


Air  temperature 

TV,  = 

10°  C 

Sea-surface  temperature 

T*  = 

o°c 

Relative  wind  speed 

V  = 

40  kt 

Wave  height 

= 

3  m 

Salinity 

s 

32V,. 

Relative  humidity 

RH  = 

90V 

Figure  48  indicates  that  the  most  significant  variables  in  determining  the  icing  severity  are 
the  air  temperature  and  the  relative  wind  speed.  The  sea-surface  temperature,  the  wave  height  and 
relative  humidity  are  seen  to  lie  of  more  moderate  significance,  while  the  salinity  of  the  sea  water  is 
shown  to  have  the  least  effect  on  rate  of  icing  of  these  six  variables. 

The  results  demonstrated  in  Figure  48  therefore  suggest  that,  in  view  of  the  quality  of  the 
available  data,  a  simplification  may  be  made  to  the  Rate  of  Icing  equation  by  assuming  a  constant 
salinity  and  a  constant  wave  height  (i.e.  a  constant  water  content,  say  0.0005  kg/nv’ ).  Further  simplifica¬ 
tion  may  be  made  by  assuming  that  the  icing  surface  temperature  is  identical  to  the  freezing  tempera¬ 
ture  of  the  sea  water.  If  then  the  icing  surface  temperature  is  taken  as  1.8°t\  there  results  the  sim¬ 
plified  equation: 


Ni*  =  0.013  (  t(1  -  1.8)V  +  0.0369V11  8  { ( - 1;1  1.8)  +  17.3(0.5354  0.9  e.,)' 


in  which  only  Td  has  to  be  determined  from  Equation  (B.14)  and  ea  by  Equation  (B.28). 

The  application  of  this  equation  to  the  data  of  Table  VIII  results  in  a  plot  not  significantly 
different  from  Figure  47  (i.e.  a  correlation  coefficient  of  0.57  and  the  best  fit  line  through  the  origin 
having  a  slope  of  0.74). 


9.0  CLOSURE 

This  report  has  collated  and  reviewed  the  work  in  the  programme'  on  the  icing  of  fishing 
trawlers  that  the  Low  Temperature  Laboratory  has  done  over  the  eleven-year  period  1969  to  1979. 
Because  of  budgetary  restraints  this  work  has  boon  very  low  key  in  nature. 
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Initially,  some  simple  icing  tests  (Ref.  17)  on  cylinders  of  various  diameters  were  made  in  the 
icing  wind  tunnel  to  demonstrate  the  effects  of  air  temperature  and  cylinder  size  on  the  resulting  ice 
formations.  The  main  conclusion  from  these  tests  was  that  the  icing  efficiency  of  a  cylinder  increases 
both  with  decreasing  diameter  and  with  decreasing  temperature,  a  result  that  demonstrates  the  need  to 
reduce  small  scale  components  such  as  spars,  rigging,  etc.  to  a  minimum,  and  where  they  cannot  he 
eliminated  altogether,  wherever  possible  to  let  one  larger  diameter  component  serve  the  purpose  of 
several  smaller  ones. 

Among  a  number  of  simple  methods  tested  as  possible  means  of  reducing  the  icing  hazard 
(Ref.  27),  no  passive  method  was  found  to  significantly  ease  the  removal  of  ice  from  a  surface.  However, 
an  inflatable  rubber  de-icing  blanket  was  found  to  be  remarkably  effective  in  removing  ice  from  both 
cylindrical  and  flat  surfaces. 

Icing  report  forms  have  been  distributed  since  1969  to  vessels  operating  out  of  eastern  Cana¬ 
dian  ports.  The  data  obtained  by  this  means  has  helped  to  establish  the  type  of  weather  systems 
responsible  for  icing  conditions,  the  geographical  extent  of  the  occurrence  of  icing,  and  statistics  on 
the  severity  of  icing  encounters  in  the  study  area. 

The  geographical  extent  varies  from  month  to  month,  Hut  at  its  most  extensive  covers  the 
whole  of  the  continental  shelf  off  Eastern  Canada  as  well  as  the  Labrador  Sea.  The  type  of  weather 
system  mostly  responsible  for  severe  icing  in  these  areas  is  a  well-developed  low  pressure  system,  in  the 
wake  of  which  cold  Arctic  air  is  brought  down  over  the  fishing  grounds  by  vigorous  northerly  or  north¬ 
westerly  winds. 

The  reports  show  that  during  January  and  February  icing  conditions  of  light  or  greater 
intensity  occur  on  the  average  one  day  out  of  two  somewhere  in  the  study  area,  while  in  bad  years  this 
may  increase  to  nearly  three  days  in  four.  In  terms  of  the  amount  of  ice  accumulated,  it  is  shown  that 
on  approximately  8  percent  of  voyages  during  the  winter  season,  ice  accumulations  in  excess  of  the 
allowance  recommended  by  IMCO  in  calculating  the  vessel's  stability  were  experienced. 

Finally,  an  analytical  approach  to  the  prediction  of  the  rate  of  ice  build-up  on  a  fishing 
trawler  is  presented.  This  analysis  (with  possible  minor  modifications)  has  application  also  to  the  icing 
of  marine  structures. 
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TABLE  I 


(Reproduced  from  Ref.  38) 


Beaufort 
Force  of 
Wind 


Theoretical 
Max.  Wave 
Height  (ft), 
Unlimited 
Duration 
and  Fetch 


Duration  of  Winds  (hours) 
with  Unlimited  Fetch,  to 
Produce  Percent  of  Max. 
Wave  Height  Indicated 


Fetch  (n.  miles),  with 
Unlimited  Duration  of 
Blow,  to  Produce  Percent 
of  Max.  Wave  Height 
Indicated 


50'/, 

75  % 

90'/ 

50'/, 

759) 

90'/, 

3 

2 

1.5 

5 

8 

3 

13 

25 

5 

8 

3.5 

8 

12 

10 

30 

60 

7 

20 

5.5 

12 

21 

22 

75 

1 50 

9 

40 

7 

16 

25 

55 

150 

280 

11 

70 

9 

19 

32 

85 

200 

450 

u<iCkUU*i  PAufc  B1AMK-M0T  FI 


4 
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TABLE  II 

ICING  QUESTIONNAIRE  RESPONSE  BY  YEAR 


Number  of  Vessels  Responding  and  Number  of  Reports 
(including  Null  Reports)  Received^ 


Year 

Fishing 

Vessels 

Government 

Vessels’ 

Freighters, 
Tankers,  etc. 

Ferry 

Services 

A 

B 

A 

B 

A 

B 

A  B 

1968-69 

2 

2 

3 

3 

1 

1 

- 

1969-70 

13 

28 

- 

- 

1 

1 

- 

1970-71 

29 

89 

- 

- 

- 

- 

1971-71 

24 

89 

1 

1 

3 

13 

1  3 

1972-732 

- 

- 

- 

-- 

1  5 

1973-74 

29 

119 

-- 

- 

1 

2 

1  1 

1974-75-’ 

5 

8 

15 

43 

8 

8 

1  40 

1975-76 

38 

172 

9 

29 

1 

2 

2  52 

1976-77 

35 

124 

10 

28 

3 

5 

4  70 

1977-78 

48 

182 

17 

76 

1 

2 

1  4 

1978-79 

22 

73 

13 

20 

Totals 

886 

200 

34 

175 

Notes: 

1  Government  vessels  include  ice  breakers,  survey  vessels,  supply  vessels,  fishery  patrol  vessels,  et< 

2  In  year  1972  73  reports  from  Newfoundland  were  lost. 

■  I  In  year  197  1  75  Newfoundland  fishing  fleets  were  on  strike 


TABLE  III 


ICING  SEVERITY  STATISTICS  FOR  THE  WHOLE  STUDY  AREA 


No.  of  Days  Icing  Severity  Reported  as: 


Light 
or  Greater 

Moderate 
or  Greater 

Heavy 

December 

7-year  total 

32 

26 

0 

worst  year 

10 

10 

3 

yearly  average 

1.6 

3.7 

1.3 

standard  deviation 

3.5 

3.2 

1.0 

daily  expectation 

1  1.7'; 

12.0'; 

1.15 

January 

7 -year  total 

OH 

70 

21 

worst  year 

23 

15 

6 

yearly  average 

1  1.0 

10.0 

3.0 

standard  deviation 

1.5 

3.0 

1.6 

daily  expectation 

15.2'; 

32.2'; 

0.75 

February 

7 -year  total 

101 

70 

32 

worst  year 

20 

If) 

7 

yearly  average 

1  1.1 

10.0 

1.6 

standard  deviation 

3.1 

3.7 

2.1 

daily  expectation 

51.0'; 

35  r; 

16.25 

March 

7 -year  total 

71 

17 

11 

worst  year 

27 

17 

;  > 

yearly  average 

10.1 

7 

2.0 

standard  deviation 

HO 

5.8 

1.0 

daily  expectation 

32.7'; 

21 .75 

6.55 

*  + 

*  *  * 

*  *  * 

1  December  —  .'(1  March 

7  year  total 

302 

2 1  3 

76 

worst  year 

72 

50 

15 

yearly  average 

13.1 

30  1 

10.0 

standard  deviation 

13  0 

10  5 

3.  1 

daily  expectation 

35  65 

25  r; 

0.05 
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TABLE  IV 

ICING  SEVERITY  STATISTICS  FOR  1CNAF  DIVISION  3L 


No.  of  Days  Icing  Severity  Reported  as: 


Light 
or  Greater 

Moderate 
or  Greater 

Heavy 

December 

7 -year  total 

7 

5 

5 

worst  year 

3 

3 

3 

yearly  average 

1.0 

0.7 

0.7 

standard  deviation 

1.1 

1.2 

1.2 

daily  expectation 

3.2'r 

2.3'; 

2.3' ; 

January 

7  -year  total 

33 

24 

9 

worst  year 

11 

8 

3 

yearly  average 

1.7 

3.4 

1.3 

standard  deviation 

2.9 

2.1 

1.3 

daily  expectation 

15. 2'; 

li. r; 

i.r; 

February 

7  -year  total 

30 

li 

2 

worst  year 

H 

i 

i 

yearly  average 

1.3 

1.6 

0.3 

standard  deviation 

1.9 

1.4 

0.5 

daily  expectation 

15.2'; 

5.6'; 

i  .0'; 

March 

7-year  total 

27 

21 

-1 

worst  year 

12 

8 

2 

yearly  average 

3.9 

3.0 

0.6 

standard  deviation 

3.9 

2.9 

0.7 

daily  expectation 

12.4'r 

9.7r; 

1.8''; 

*  *  * 

*  * 

*  *  * 

1  December  —  31  March 

7 -year  total 

97 

61 

20 

worst  year 

34 

20 

5 

yearly  average 

13.9 

8.7 

2.9 

standard  deviation 

8.8 

5.3 

1.8 

daily  expectation 

11.4'; 

7.2'; 

2,r; 
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TABLE  V 

ICING  SEVERITY  STATISTICS  FOR  ICNAF  DIVISION  3P 


No.  of  Days  Icing  Severity  Reported 

as: 

Light 

or  Greater 

Moderate 
or  Greater 

Heavy 

December 

7-year  total 

4 

3 

0 

worst  year 

2 

2 

0 

yearly  average 

0.6 

0.4 

0 

standard  deviation 

0.7 

0.7 

0 

daily  expectation 

1.8',/ 

1.4'; 

O'; 

January 

7 -year  total 

53 

38 

10 

worst  year 

n 

8 

3 

yearly  average 

7.6 

5.4 

1.4 

standard  deviation 

3.0 

2.7 

1.2 

dialy  expectation 

24.4'; 

17.5'; 

4.6'; 

February 

7-year  total 

59 

41 

18 

worst  year 

16 

12 

4 

yearly  average 

8,1 

5.9 

2.6 

standard  deviation 

3.6 

3.1 

1.3 

daily  expectation 

29. H'; 

20.7'; 

9.r; 

March 

7-year  total 

38 

23 

i 

worst  year 

14 

7 

3 

yearly  average 

5,1 

3.3 

0.6 

standard  deviation 

5.2 

2.8 

1.0 

daily  expectation 

1 7.5'; 

10.6'; 

1  ,8'V 

*  * 

*  *  * 

*  *  * 

1  December  —  31  March 

7-year  total 

1 54 

105 

32 

worst  year 

37 

27 

7 

yearly  average 

22.0 

15.0 

4.5 

standard  deviation 

8.3 

6.3 

1.7 

daily  expectation 

18  r; 

12.4'-; 

3.8'; 

SUMMARY  OF  ICING  ENCOUNTERS  IN  TERMS  OF  IMCO  STANDARD  ICING  ALLOWANCE 
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TABLE  VII 

NUMBER  (AND  PERCENTAGE)  OF  YEARS  FOR  WHICH  AT  LEAST  ONE  OCCURRENCE  OF 
ICING  OF  SEVERITY  GREATER  THAN  SLIGHT  HAS  BEEN  REPORTED  FOR  THE 
DESIGNATED  ICNAF  DIVISION  IN  EACH  OF  THE  WINTER  MONTHS 


ICNAF 

Division 

Nov. 

Dec. 

Jan. 

Feb. 

Mar. 

Apr. 

IF 

1  (13%) 

2G 

1  (13%) 

2H 

1  (13%) 

2  (25%) 

1  (13%) 

2J 

3  (38%) 

3  (38%) 

3  (38%) 

1  (13%) 

3K 

2  (25%) 

1  (13%) 

4  (50%) 

4  (50%) 

2  (25%) 

3L 

1  (13%) 

2  (25%) 

7  (88%) 

7  (88%) 

4  (50%) 

3M 

3N 

1  (13%) 

2  (25%) 

30 

1  (13%) 

2  (25%) 

5  (63%) 

2  (25%) 

3P 

3  (38%) 

8  (100%) 

8 (100%) 

6  (75%) 

2  (25%) 

4R 

1  (13%) 

2  (25%) 

5  (63%) 

4  (75%) 

1  (13%) 

2  (25%) 

4S 

3  (38%) 

4  (50%) 

1  (13%) 

1  (13%) 

4T 

2  (25%) 

2  (25%) 

3  (38%) 

1  (13%) 

4Vn* 

4  (50%) 

8  (100%) 

6  (75%) 

3  (38%) 

2  (25%) 

4Vs* 

4W 

4  (50%) 

4  (50%) 

7  (88%) 

2  (25%) 

4X 

3  (38%) 

5  (63%) 

4  (50%) 

3  (38%) 

5Y 

1  (13%) 

1  (13%) 

2  (25%) 

1  (13%) 

Based  on  8  years  of  data 


* 


The  boundary  between  4Vn  and  4Vs  is  here  taken  to  be  the  parallel  44°  10'  N  instead  of  45°  40’  N. 
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TABLE  VIII 

DATA  USABLE  FOR  ICING  RATE  STUDIES 
FROM  SEASONS  1977-78  AND  1978-79 


(1) 

Data 

Set 

No. 

(2) 

Vessel 

Length 

m 

(3) 

Rel.  Wind 
Velocity 
kt 

(4) 

Air 

Temp. 

°C 

(5) 

Sea-Surf. 

Temp. 

°C 

(6) 

Wave 

Height 

m 

(7) 

Sea 

Salinity 

%o 

(8) 

Icing 

Rate 

mm/h 

(9) 

Calculated 
Icing  Rate 
(Eq.  B.3/c) 
mm/hr 

1 

51.5 

55 

-6 

+1 

5* 

32.25 

8.5 

5.1 

2 

51 

35 

-3.5+ 

-0.5 

1.2 

31 

1.6 

2.1 

3 

51 

35 

-6  + 

+4 

3* 

32.35 

0.79 

3.3 

4 

51 

48 

-3 

-0.5 

4.5* 

34 

4.2 

1.6 

5 

51 

32 

-5 

-1 

1 

33.75 

1.3 

3.0 

6 

51 

56 

-4 

-2 

3* 

33.25 

0.75 

4.2 

7 

51 

36 

-7  + 

+1 

3* 

32.25 

10.9 

5.3 

8 

51 

39 

-4 

0 

3* 

32.7 

3.2 

2.6 

9 

51 

32 

-10  + 

+2 

9 

33.75 

4.8 

7.3 

10 

47 

36 

-2 

0 

0.6* 

31.5 

4.8 

0.65 

11 

47 

25 

-5 

0 

2.4* 

34 

0.85 

2.7 

12 

47 

40 

-9  + 

0 

2.4 

32.5 

0.71 

7.9 

13 

47 

40 

-6.5 

+0.5 

3* 

32.75 

2.8 

5.2 

14 

47 

50 

-4 

-1 

3* 

32.3 

5.4 

3.5 

15 

47 

55 

-8 

+1 

3* 

32.3 

1.6 

8.2 

16 

46 

38 

-8  + 

+1.5 

2 

31 

3.8 

6.3 

17 

44 

45 

-14 

-0.5 

3* 

32.75 

12.7 

14.1 

18 

44 

35 

-5.5 

+4 

1 

32.1 

3.0 

3.3 

19 

44 

28 

-3 

+1 

2* 

32 

1.8 

1.2 

20 

44 

28 

-10 

0 

2* 

31 

7.3 

6.8 

21 

42 

37 

-10 

-0.5 

2* 

33.5 

3.6 

8.2 

22 

42 

26 

-5 

0 

1* 

31 

2.6 

2.6 

23 

42 

25 

-5 

0 

1* 

33.5 

2.1 

2.4 

24 

42 

25 

-5 

+3 

3* 

34 

1.9 

2.2 

25 

42 

25 

-6 

0 

1.5 

31 

8.0 

3.4 

26 

42 

33 

-10 

+1.5 

2 

32.5 

3.5 

7.3 

27 

42 

38 

-12 

+2 

1.5 

32.25 

6.5 

9.5 

28 

42 

30 

-10 

0 

1.2 

32 

6.9 

6.6 

29 

42 

39 

-16 

0 

1.2 

33.75 

10.6 

12.8 

30 

38 

41 

-7 

-1 

3* 

32.25 

1.6 

6.4 

31 

34 

26 

-10  + 

0 

3* 

31 

2.4 

6.8 

32 

34 

32 

-7  + 

-1 

3* 

32 

2.2 

5.4 
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TABLE  VIII  (Cont’d) 

DATA  USABLE  FOR  ICING  RATE  STUDIES 
FROM  SEASONS  1977-78  AND  1978-79 


(1) 

Data 

Set 

No. 

(2) 

Vessel 

Length 

m 

(3) 

Rel.  Wind 
Velocity 
kt 

(4) 

Air 

Temp. 

°C 

(5) 

Sea-Surf. 

Temp. 

°C 

(G) 

Wave 

Height 

m 

(7) 

Sea 

Salinity 

%  o 

(8) 

Icing 

Rate 

mm/h 

(9) 

Calculated 
Icing  Rate 
(Eq.B.3/c) 
mm/hr 

33 

51 

46 

-6  t 

0 

2 

32.25 

4.0 

5.3 

34 

51 

21 

-10 

-0.5 

3* 

33.25 

3.9 

5.8 

35 

47 

29 

-5 

+0.5 

2 

33.5 

2.8 

2.9 

36 

44 

31 

-6 

0 

2* 

32.3 

1.9 

4.0 

37 

42 

35 

-6 

-0.5 

3* 

32.5 

2.1 

4.6 

38 

42 

19 

-4 

0 

2* 

32.2 

1.6 

1.7 

39 

38 

31 

-5 

+2 

3 

31 

2.5 

2.9 

Data  in  Columns  (3),  (4), (6)  and  (8)  derived  from  Ship  Icing  Reports 

t  These  air  temperatures  inferred  from  other  vessels  in  same  general  vicinity. 

*  These  wave  heights  inferred  from  sea  state  description. 


IN  DATA  SETS  1  TO  32  THE  RELATIVE  WIND  DIRECTION  IS  0°  ±45°. 
DATA  SETS  33  to  39  HAVE  RELATIVE  WIND  DIRECTION  >  45°. 


DISTANCE  BETWEEN  UPPER  AND  LOWER  TANGENT  TRAJECTORIES  IN 


>rtECE£>iNtt  PAGE  BUNK -WOT  F1u>&) 


FIG.  1 :  WATER  DROP  TRAJECTORIES  IIM  THE  VICINITY  OF  A  CIRCULAR  CYLINDER 


FIG.  2  ICE  ACCRETIONS  ON  15.2  CM  HORIZONTAL  CYLINDER  AFTER  1  HOUR  EXPOSIJRI 


CCRETIONS  ON  15  2  CM  VERTICAL  CYLINDER  AFTER  1  HOUR  EXPOSURE 

RUN  NO  16 
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~r 


E 

o 

cn 


10-20%  LESS 
THAN  MIDSECTION 


MIDSECTION 


10-20%  GREATER 
THAN  MIDSECTION 


FIG.  4:  DIMENSIONS  OF  ICE  ON  VERTICAL  AND  HORIZONTAL  CYLINDERS 


ICING  EFFICIENCY  % 
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AIR  TEMPERATURE  °C 


FIG.  5:  EFFECT  OF  TEMPERATURE  ON  ICING  EFFICIENCY 


FIG.  8.  TWO  PHASE  THERMOSYPHON 


G  A 1  j  ij  £ 

* —  2  -  A  A  Y  0  P  f  R 

^--PRESSURE  PFl  ff  v A i  v  t 
PRESSURE  REDUCING  VALVE 

U>!  CONTROLGEAR 


FIG.  10:  DETAILS  OF  PNEUMATIC  DE  ICER 


TEMPERATURE:  IOC  ICE  THICKNESS:  UP  TO  9  cm  DE-ICER  AIR  PRESSURE:  70  kPa 


r 


(a)  AT  CONCLUSION  OF  ICING  (b)  AFTER  TWO  INFLATIONS  OF  DE' ICE R 


(c)  AFTER  TWO  INFLATIONS  OF  DE  ICER 


Id)  AFTER  THIRD  INFLATION 


FIG.  14:  OUTDOOR  TEST  NO.  2 


ICING  TEMPERATURE:  15°C 
DE  ICED  AT  -19°C  AFTER  17  HOURS  COLD  SOAK 


(a)  AT  CONCLUSION  OF  ICING 


(bl  AFTER  OPERATION  OF  DE  ICERS 


FIG.  15:  OUTDOOR  TEST  NO.  6 


ICING  TEMPERATURE  15°C 
DEICED  AT  12°C 


FIG.  16:  PRINCIPLE  TRACKS  OF  INTENSE  STORMS  FOR  JANUARY, 
FEBRUARY  AND  MARCH  1963-67 


(FROM  ARCHIBALD,  Ref.  36) 


FIG.  17.  MEAN  SEA-SURFACE  TEMPERATURE  (  F)  FOR  FEBRUARY 

(FROM  OCEANOGRAPHIC  ATLAS  REF  371 
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FIG.  18;  FREEZING  TEMPERATURE  OF  SEA  WATER  AS  A  FUNCTION  OF  SALINITY 
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SHIP  ICING  RETORT 


HI  l  AM  HtPOHI  UN  IMIS  I  OHM  ANY  INSTANCE  Of-  l<  t  I  (JRM  A 1 1  UN  UN  VISbtl  DUHINC.  VOVA!.*,  H 
MOW  1  MAN  ONI  OO  UHt’lNO  O*  lUNG  ON  VOYAOF.  Ubl  btPARAH  I  OHM  I  QH  IAC.H  INC  UONUH 
IT  I-  IMPORTANT  1()  <  OMPUU  A  I  OHM  MJH  t  AC  H  VOYAcU  IVFN  II  NO  IC_:ifS|( »  QC.CURb. 


f  OMPl  t  Tl  [I  »  ORM‘.  SMOUI  0  R{  MAH  1 1)  TO  .  =  >*  » » ma  >•  hi  ,  ah.  .ha  »i  .h>  i ..  .-mi  in  -  .i  mu  ••  a  m,  a  UpAIiiov 

•lA'il'jA,  HI  M  AMi  Mill  N-  •  •  '  t  AAA  •  1%  1  AH  ,  .1  *  <  A  Hi 


FIG.  19:  SHIP  ICIN?  REPORT  FORM 


FIG.  22:  DESIGNATED  SEA  AREAS  OF  THE  INTERNATIONAL  COMMISSION  FOR 
THE  NORTHWEST  ATLANTIC  FISHERIES  (ICNAF) 


FIG.  23:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  NOVEMBER 


NUMBER  OF  YEARS  OF  DATA  8 


FIG.  24:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  DECEMBER 

NlJMIjfH  Of"  Vf  ARSOF  DA1A  8 


■  ()1 


FIG.  25:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  JANUARY 


NUMBS  R  OF  YF  AKS  OF  DA  I  A  8 


FIG.  26:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  FEBRUARY 


NUMBER  OF  YEARSOF  DATA  8 


FIG.  27:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  MARCH 

NUMBER  OF  YEARS  OF  DATA:  8 


FIG.  28:  PERCENTAGE  OF  YEARS  FOR  WHICH  ICING  OF  SEVERITY  GREATER  THAN 
SLIGHT  HAS  BEEN  REPORTED  IN  THE  VARIOUS  ICNAF  DIVISIONS  FOR  THE 

MONTH  OF  APRIL 

NUMBER  OF  YEARS  OF  DATA:  8 


-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  29:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°  F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  O'C  OR  LESS 
FOR  THE  MONTH  OF  NOVEMBER 


A 
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-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

- 10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  30:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  DECEMBER 


+4.4  C  MEAN  SEA  SURFACE  ISOTHERM 


FIG.  31 


L 


10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 
APPROXIMATE  MEAN  ICE  LIMIT 


:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°  F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  JANUARY 


-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- APPROXIMATE  MEAN  ICE  LIMIT 


FIG.  32:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  O  C  OR  LESS 
FOR  THE  MONTH  OF  FEBRUARY 


r 


-  +4.4°C  MEAN  SEA  SURFACE  ISOTHERM 

-  10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 

- - APPROXIMATE  MEAN  ICE  LIMIT 

FIG.  33:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  MARCH 


r 


+4.4  C  MEAN  SEA  SURFACE  ISOTHERM 


10%  PROBABILITY  OF  AIR  TEMPERATURE  0°C  OR  LESS 
APPROXIMATE  MEAN  ICE  LIMIT 


FIG.  34:  MEAN  SEA  SURFACE  TEMPERATURE  OF  40°F  (4.4°C)  and  10%  FREQUENCY 
OF  OCCURRENCE  OF  SURFACE  AIR  TEMPERATURE  OF  0°C  OR  LESS 
FOR  THE  MONTH  OF  APRIL 


X  -  VESSELS  EXPERIENCING  ICING 
■  VESSEL  NOT  EXPERIENCING  ICING 

FIG.  35:  SURFACE  CHART  SHOWING  METEOROLOGICAL  SYSTEM  RESPONSIBLE  FOR 
ICING  OCCURRENCES  ON  14  FEBRUARY,  1970 


X  -  VESSELS  EXPERIENCING  ICING 
•  -  VESSEL  NOT  EXPERIENCING  ICING 

FIG.  36:  SURFACE  CHART  SHOWING  METEOROLOGICAL  SYSTEM  RESPONSIBLE  FOR 
ICING  OCCURRENCES  ON  17  FEBRUARY,  1970 


X  VESSELS  EXPERIENCING  ICING 

■  VESSELS  NOT  EXPERIE  NCING  ICING 


SHELTERING  IN  LEE  OE  LAND 

FIG.  37  SURFACE  WEATHER  CHART  FOR  1200  GMT,  27  JANUARY,  1972 


k 


X  VESSELS  EXPERIENCING  ICING 


■  VESSEL  IN  PACK  ICE  -  NOT  E XPE  RIE  NCING  ICING 

FIG.  38:  SURFACE  WEATHER  CHART  FOR  0600  GMT,  12  MARCH,  1974 


X  -  VESSELS  EXPERIENCING  ICING 


FIG.  39:  SURFACE  WEATHER  CHART  FOR  0000  GMT,  14  FEBRUARY,  1979 


X  VESSELS  EXPERIENCING  ICING 


FIG.  40:  SURFACE  WEATHER  CHART  FOR  1800  GMT,  15  FEBRUARY,  1979 


(a)  ICE  BUILD-UP  ON  RAILS,  ANLhOR  AND  FORESTAY 


(b>  ICE  BUILD  UP  ON  WHEELHOUSE  FRONT 


FIG.  41:  ICING  OF  STERN  TRAWLER  26  27  JANUARY  1972 


relative  wind  speed  -  KNOTS 


1 

-  79  -  ] 


AIR  temperature  -  °c 


FIG.  43:  EFFECT  OF  AIR  TEMPERATURE  AND  RELATIVE  WIND  SPEED  ON 
REPORTED  ICING  SEVERITY  FOR  FISHING  VESSELS  STEAMING 
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FIG.  44:  EFFECT  OF  AIR  TEMPERATURE  AND  RELATIVE  WIND  SPEED  ON 
REPORTED  ICING  SEVERITY  FOR  FISHING  VESSELS  AT  LOW  SPEEDS 

(4  KNOTS  AND  BELOW) 


AIR  TEMPERATURE 


BEAUFORT  FORCE  8 


BEAUFORT  FORCE  6-7 


BEAUFORT  FORCE  9-10  BEAUFORT  FORCE  11-12 


SEA -WATER  TEMPERATURE  (°C) 


DEGREE  OF  ICING  0 
1 
2 

3 

4 


NONE 
LIGHT 
MODERATE 
SEVERE 
VERY  SEVERE 


1  3  CM  ?A  HR 
A  6  CM  24  HR 
/  14  CM  24  HR 

•  1  E>  CM  '  24  HR 


FIG.  45:  MERTINS'  CHARTS  OF  ICING  ON  FISHING  VESSELS  AT  SLOW  SPEEDS 
IN  WINDS  OF  BEAUFORT  FORCE  6  12 

(REPRODUCED  1  ROM  HI  E  4b) 


0  2  4  6  8  10  12  14 

COMPUTED  RATE  OF  ICING  -  mm/h 


FIG.  47:  COMPARISON  BETWEEN  CALCULATED  RATE  OF  ICING  AND  THAT  ESTIMATED 

FROM  TRAWLER  ICING  REPORTS 


BROKEN  LINE  REPRESENTS  BFST  LEAST  SQUARES  FIT  THROUGH  THE  ORIGIN.  SLOPE  t)  8? 


ICING  , 
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SIMPLI I'lKI)  KORMM.AK  K)H  TOT Al,  COLLECTION  I  I  I  ICII  NOV 
()!•  CYI.INDKRS  AND  RI.CTAMU  I.AK  DODDS 


In  an  attempt  to  simplify  the  calculation  of  the  droplet  collection  efficiency  of  cylinders  and 
rectangular  bodies  and  to  obviate  the  use  of  curves  such  as  presented  in  References  14  and  19  the  fol¬ 
lowing  formulae  are  proposed.  These  formulae  have  been  derived  for  the-  rather  limited  range  of 
atmospheric  conditions  appropriate  to  trawler  icing  and  the  icing  of  static  structures. 

The  properties  of  air  at  -  5°C  and  100  kl’a  (750  mm  1  Ig)  have  been  used,  and  the  deviations 
from  these  conditions  that  are  to  be  expected  in  this  context  have  negligible  effect  on  the  final  result. 
Thus  only  the  variables  of  velocity,  droplet  diameter,  and  cylinder  diameter  (or  body  width )  arc  in¬ 
volved  in  the  formulae  and  are  combined  into  a  single  parameter,  if  (zeta),  of  the  form: 

f  =  V1'  d‘<  IT  (All 


when*  V  =  velocity 

d  =  droplet  diameter 

I)  =  cylinder  diameter  (or  body  width) 

The  indices  p,  q  and  r  must  be  chosen  such  that  the  families  of  curves,  resulting  from  the  use 
of  the  K  and  (or  Reo  )  parameters  as  used  in  References  14  and  49  all  collapse  to  a  single  curve  when 
the  collection  efficiency,  K  ,  is  plotted  against  f. 

It  was  found  that  this  requirement  was  satisfied  over  a  representative  range  of  values  of  V,  d 
and  I)  (with  the  exception  of  the  simultaneous  combination  of  very  small  values  of  all  three  parameters) 
by  the  following  values  of  the  indices: 

p  =  ().(> 

q  =  1.(5 
r  =  1 

V"  (,  ,,1.0 

whence  f  I  A. la) 

I) 


I 


The  units  used  throughout  in  evaluation  f  are: 


for  V  m/s,  for  d  pm,  and  for  I)  m 


I  o  conform  ’o  the  obvious  requirements  that  when  V  or  d  is  very  large,  or  when 
small  I  i.e.  when  f  •  '■«),  the  collection  efficiency.  Km  ,  approaches  unit  v.  and  that  K 
<T  •  0.  a  relat  ion  ship  hot  wren  f  and  Km  of  the  following  form  was  adopted 


I )  is  very 
0  when 


-  8fi  - 


K,„  -  -  for  f  J 

cv 

1  +  ~  - - 

tf  fi)" 


K 


m 


0 


fur  ?  •  p 


(A  2t 


where  a,  p  and  n  are  constants. 

Cylinder  Collection  Efficiency 

Evaluation  of  a,  /j  and  n  for  the  cylinder  collection  efficiency  data  presented  in  Reference  1  1 
leads  to  the  following  relations: 


1  + 


1 

3200 
?  3200 


f  3200 
f  +  27000  ' 


3200 


0 
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This  curve  is  presented  in  Figure  A.l,  together  with  actual  values  of  Kr  evaluated  using 
Reference  11  over  the  following  range  of  values: 

V  from  3  m/s  to  (50  m/s 
d  from  20  pm  to  1000  pm 
I)  from  0.03  m  to  1  m 

As  the  points  indicate,  the  values  of  E|n  are  collapsed  w el]  by  Equat ion  (  A.  1  a |  onto  a  single 
curve,  with  the  exception  of  a  few  points  which  lie  below  the  curve,  part  ieularly  at  its  lower  end.  and 
which  correspond  to  combinations  of  the  lower  values  of  the  parameters  Y,  d  and  I ). 

The  curve  given  by  Equation  (  A. 3)  is  seen  to  correspond  to  better  than  0.0  1  in  E(ll  value  to 
the  main  group  of  points,  while  the  cx'raneous  points  hi'  no  more  than  0.08  m  E  value  from  the 
curve. 


Rectangular  Body  Collection  Efficiency 


1'he  rectangular  body  collection  efficiency  data  of  Reference  10  were  used,  anil  led  to  the 
following  relations: 
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Figure  A.l  presents  tills  curve  toget  her  with  values  of  E.  e\  a  licit  ed  Using  R  cl  ere  III  I  10  over 
the  same  range  ,»f  the  variables  V  an  d  as  for  the  cylinder  data,  however  the  bo,  lv  w  id  lit  1 1 ) )  was 
extended  to  3 m  Die  R,„  data  collapses  onto  a  single  curve  rat  her  better  than  for  the  cylinder  case, 
except  win  n  the  droplet  Reynolds  Number  i  Re  |  is  less  than  about  20.  in  win,  h  case  the  points  fall 
below  tin-  mam  grouping  of  poults.  I  In  -  <  urve  i  Eq.  i  A  In  agrees  well  wnii  the  main  ln„|\  of  point  - 


COLLECTION  EFFICIENCY  OF  CYLINDER  AND  RECTANGULAR  BODY 
AS  A  FUNCTION  OF  PARAMETER  ? 
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APPENDIX  B 

DERIVATION  OF  EXPRESSION  FOR  PREDICTING 
THE  ICING  R  ATE  OK  FISHING  TRAW  LERS 

B.1.0  HEAT  BALANCE  AT  THE  ICING  SURFACE 

Taking  only  the  primary  heat  transfer  processes  acting  at  the  icing  surface  into  consideration, 
and  assuming  that  the  ice  formation  is  a  continuous  steady-state  process,  an  equilibrium  heat  balance 
at  the  icing  surface  may  be  formulated  as  follows: 

q,  +  qw  +  qt.  +  qt,  +  qa  =  0  (B.ll 

where  qt  =  latent  heat  of  freezing  a  certain  fraction  of  the  impinging  water. 

qtt  =  heating  (or  cooling)  of  the  impinging  water  to  the  equilibrium  surface 

temperature, 

q  =  heat  transfer  by  convection  to  the  surrounding  air, 

q0  =  heat  transfer  by  evaporation  to  the  surrounding  air. 

q.,  =  heat  conducted  to  or  from  the  icing  surface  through  the  underlying 

structure. 

Each  of  these  terms  will  initially  be  considered  individually  in  the  following  sub-sections. 
B.1.1  The  Latent  Heat  Term,  q, 

The  rate  of  water  impingement  on  unit  area  of  a  surface  may  be  expressed  as: 

Rw  =  KwV  (  B.2 1 

where  K  is  the  collection  efficiency. 

w  is  the  liquid  water  content  of  the  air. 

V  is  the  velocity  of  the  air  relative  to  the  object  or  surface. 

I'he  collection  efficiency  expresses  the  degree  to  which  the  water  drop.'  m  the  air  impinge  on 
an  object  in  I  heir  path .  and  are  not  deflected  bv  the  air  How  around  the  object .  I'he  collect  ion  el  fan  cm  a 
increases  with  relative  velocity,  and  with  drop  size,  and  decreases  with  the  size  of  tin-  object  on  w  Inch 
tile  drops  impinge.  Because  the  Size  of  till'  water  drops  involved  111  slco  icing  due  to  sea  spray  Is  large 
I  1  nun  ).  they  will  he  deflected  little  and  t  he  collect  ion  efficiency  will  be  assumed  in  be  IDO 

llte  distribution  of  drop  sizes  m  the  spray  is  unknown  and  vanes  with  wind  speed,  -h  |  ■  -peed 
t  V  pe  and  Size  of  vessel .  the  local  loll  on  the  Vessel,  and  with  height  above  the  sea  surface.  Data  quoted 
bv  Borisenkov  i  Ref.  M  i  suggests  the  drop  diamei  er  varies  f  r<  >m  1  .0  to  :!.f>  mm  (the  average  is  \> .  I  mini  a,, 
a  result  of  spray  mg  fr<  mi  swell  waves  I  low  ever,  in  the  bow  of  the  ship  the  w  a  ter  will  lie  encountered  in 
the  form  of  "giisle  or  “ease.i.ies"  for  vvhteli  i  lie  concept  of  drop  size  seem-  ;na |  | >ro print ■  uev  «*r:  He!*  - 
for  Dio  w  ant  of  a  i  id!  er  concept ,  g  will  In  employ'd  in  1 1 1  e-  anal  v  i 

I  111'  liquid  water  content  of  the  sprav  elond  |s  rqllallv  difficult  to  quanttfv  ill  .  It  I  \  U'liqHi  w  a\ 
however.  Kaelmrm  t  Rel'  1  7  I  -Ugges1  .  t  hat  wave  height  is  the  chief  fa<  lor  governing  tin  p.namitil  \ 
w  it  1 1  drop  si/e.  it  may  be  as'1  miei  I  t  bat  lb.  vv.iti  r  eonteiit  e  also  (ui  a!  ion  of  tin  mz<  .  I  n<  -  and  -peed 
of  t  lie  v  e.ssel.  and  al-o  a  to  net  ion  ot  i  eight  above  tin  in  sin  In  e.  I\a<  h  linn  propo.i  ;hat  for  '..ring 
-peed-,  ot  1 1  's  I- 1 1 ' 1 1  -  at  wl  In I  mg--  o|  li  •  ill  _  a  propori  an  ial  i  a  iai  loii-ii  ip  i  n  .  is-i  I  nied 
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w  =  10"1  llu  kg/nr* 

where  Htt  is  the  wave  height  in  metres. 

However,  in  the  icing  model  developed  here,  it  was  found  that  a  liquid  water  content  one 
sixth  of  that  given  by  expression  (B.3)  resulted  in  a  significant  improvement  in  correlation  with  observed 
results.  As  a  result,  the  expression  for  liquid  water  content 

w  =  1.7  X  10-4  Hu  kg/trv’  (B.dal 


has  been  adopted. 


If  a  fraction,  n,  of  the  impinging  water  freezes  on  the  surface,  then  the  rate  of  ice  growth  per 
unit  area  is: 


R,  =  nR„  =  nKwV 

and  if  is  the  latent  heat  of  freezing,  then  the  heat  released  is: 

9 1  =  <qR, 


(  11.  1  I 


i  B.r>  i 


A.  1 .2  'Hie  Sensible  Heat  Term,  qw 

If  the  equilibrium  temperature,  ts,  of  the  icing  surface  on  which  the  water  drops  impinge  is 
assumed  to  be  the  equilibrium  freezing  temperature,  t, .  of  water  of  the  appropriate  salinity  (i.e.  ignoring 
any  temperature  gradient  across  any  water  film  on  the  icing  surface),  and  the  drop  temperature  imme¬ 
diately  prior  to  impingement  is  t  then  the  heat  given  to  the  impinging  water  is: 

<1„  ~  R«  *V  'hi  -  '  J  (B.d) 

where  ew  is  tig'  specific  heat  of  the  impinging  water  >  1000 -J  kgC 

lor  s  --  35';.,  (•'  IOC. 

It  should  not  lie  overlooked  however  that  some  of  ihe  impinging  water  w  ill  likely  splash  or 
run  off  the  icing  surface  before  reaching  the  equilibrium  temperature,  tv.  Thus  a  splash  factor,  s.  should 
appropriately  be  introduced  into  Ivqtiat  ion  |  B.(> )  as  follows: 

(),,  I  I  s)  Kw  ctt  (t(|  t  J  I  li.tia ) 

However,  t  here  seems  lit  t  le  way  of  knowing  the  magi  lit  tide  of  the  factor,  s,  but  .i;iclwri'' 

I  Ref.  17  I  suggests  that  about  0.75  of  the  water  film  runs  off  without  treezmg  and  assumes  there  is  no 
contribution  whatever  to  the  heat  balance  from  this  water.  In  fact,  it  is  likely  that  tile  temperature  of 
the  water  film  is  close  to  that  of  the  ice  surface,  and  that  it  is  only  the  water  that  splashes  and  is  blow  n 
away  at  the  moment  of  impingement  that  takes  no  part  in  the  thermal  process. 


Because  of  these  uncertainties  the  form  of  Kquatlon  I  B.fil  will  be  employed  U,  ihe  follow  mi. 

analysis 


I  he  tempera!  lire  of  the  mu miging  water  drops,  t  e  a  fund  ion  of  man \  van.  hies  melu,  lim. 
the  sea -surface  I  dupe  rat  lire,  the  air  1  cm  pern  I  tire  the  drop  si/e  and  the  time  ot  flight . 
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An  expression  for  the  drop  temperature  may  be  derived  by  following  the  procedure  due  to 
Hardy  (Ref.  50). 


The  rate  of  change  of  temperature  of  the  drop  due  to  convection  is: 


Nu 


6  ka 
Pw  ^ 


(ta- 1) 


where  t 
t’a 

T 

Nu 

ka 

Pw 

Cw 

d 


instantaneous  drop  temperature 

air  temperature 

time 

Nusselt  No. 

thermal  conductivity  of  air 
density  of  water 
specific  heat  of  water 
drop  diameter 


(B.7) 


The  rate  of  evaporation  is  given  by: 


where  m 


D 

Pa 

Mw/Ma 
e,  ea 
P 


dm  Mw  /V  e\ 

—  -  Nu„dDp.  —  J  I 

mass 

diffusivity  of  water  vapour  in  air 
density  of  air 

ratio  of  molecular  weights  of  water  and  air  (0.622) 

vapour  pressure  at  surface  of  drop,  and  in  surrounding  air  respectively 

barometric  pressure 


But  the  temperature  change,  dt,  resulting  from  an  evaporative  mass  loss,  dm,  from  a  drop  of 
diameter  d  is: 


dt 


6 


Pw  ®w  ^ 


dm 


(B.9) 


where  Cv  =  latent  heat  of  vaporization  (2.5  X  106J/kg). 
Thus,  the  rate  of  change  of  temperature  due  to  evaporation  is: 


ggyDPa  Mw  /ea  -  e\ 

pwcwd2  \  p  J 


(B.10) 


Combining  Expressions  (B.7)  and  (B.10)  to  obtain  the  overall  rate  of  temperature  change: 
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dt  _  6  Nu  ka 
dr  pwcwd2 

6  Nu  ka 

=  p  c  d2  ^  "  L> 
where  X,  is  a  factor  to  allow  for  evaporation: 

ea  _  e 

X,  =  1  +  0.622  -  - 

PCp  ta  -  t 


(t,  -  t)  ll  + 


Mw£vDpa 

Makap 


G3I 


V  /I 


(B.ll) 


(B.lla) 


(B.12) 


k 

since  D  — - 

Pacp 


and  where  cp  =  specific  heat  of  air  (1.005  X  103  J/kgK). 


Since  Xt  is  a  weak  function  of  temperature.  Expression  (B.lla)  may  be  integrated  assuming 
X,  is  a  constant,  resulting  in: 


(ta  -  td )  6  Nu  ka 


(B.13) 


or  td  =  ta  +  (tw  -  ta)  e 


/  6  Nu  k3  \ 

V  PwCw«i-  *  / 


(B.14) 


where  td  is  the  water  drop  temperature  at  time  r. 


To  test  the  validity  of  the  assumption  that  the  factor  X,  may  be  treated  as  a  constant  in  the 
integration,  results  obtained  using  Equation  (B.14)  are  compared  in  Figure  B.l  with  the  results  of  a 
numerical  integration  of  Equation  (B.ll).  Also  shown  for  comparison  is  the  effect  of  ignoring  the 
evaporative  effect  (convection  only).  It  is  clear  that  the  assumption  of  a  constant  X,  does  not  result 
in  any  serious  errors  in  the  computation  of  drop  temperature. 

The  appropriate  time  or  distance  of  flight  of  a  drop  before  striking  some  part  of  the  ship  is 
uncertain,  as  is  also  a  representative  drop  diameter.  Clearly,  larger  drops  will  stay  airborne  for  shorter 
times  than  smaller  drops,  so  that,  together  with  the  increase  in  the  value  of  the  time  constant 

l in  Equation  (B.14)  resulting  from  a  larger  d,  larger  drops  experience  considerably  less 

cooling  than  smaller  drops,  a  circumstance  which  results  at  times  in  the  decks  and  lower  parts  of  a  vessel 
experiencing  little  ice  build-up  while  appreciable  ice  forms  on  rigging  and  superstructure. 


ffiw  cwd2 


,6  ka  Nu 


For  the  sake  of  this  analysis,  however,  an  effective  drop  diameter  of  2  mm  is  assumed.  It  is 
further  assumed  that  the  effective  relative  velocity  between  the  drop  and  its  surrounding  air  is  the  fall 
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velocity  of  the  drop,  which  for  a  2  mm  drop  is  6.5  m/s  (see  Ref.  51 ).  Thus,  using  the  expression  for 
the  convective  heat  transfer  for  a  sphere  given  by  McAdams  (Ref.  52)  i.e.  Nu  =  0.37  Re0  6,  the  Nusselt 
number  is  23  and  hence  the  time  constant  is  approximately  5/X,  seconds. 

The  flight  time  of  the  drops,  r,  was  expressed  in  Reference  42  as  x/V  (i.e.  the  distance  carried 
by  wind  of  velocity  V)  and  it  was  suggested  that  x  might  be  related  to  vessel  length  as  a  means  of  intro¬ 
ducing  vessel  size  into  the  analysis  however,  this  is  perhaps  not  the  appropriate  place  to  introduce  size 
since  the  effect  would  be  of  increasing  icing  rate  with  increasing  vessel  length,  a  situation  that  does  not 
always  accord  with  observation.  Perhaps  some  modification  to  the  liquid  water  content  expression 
(Eq.  (B.3))  would  be  the  more  appropriate  way  of  introducing  both  the  size  and  lines  of  the  vessel. 

More  knowledge  of  the  effect  of  size  and  lines  on  the  rate  of  icing  is  first  needed. 

In  this  analysis  the  flight  distance,  x,  is  assumed  constant  at  20  m,  whence  the  effective  flight 
time  of  the  drops,  t,  becomes  20/V,  a  figure  that  results  in  the  best  correlation  between  the  analytical 
and  observed  rates  of  icing. 

The  surface  temperature,  ts,  was  stated  earlier  to  be  the  equilibrium  freezing  temperature  of 
water  of  the  appropriate  salinity.  As  ice  forms  at  the  surface,  salt  is  rejected  from  the  ice  crystal  lattice 
so  formed,  and  the  water  film  on  the  surface  of  the  ice  manifests  an  enriched  salt  content.  Since  a 
definite  relationship  exists  between  the  salinity  and  the  freezing  temperature  of  the  brine,  it  follows 
that  the  surface  temperature,  ts,  must  be  that  appropriate  to  the  salinity  of  the  brine  on  the  surface  of 
the  ice,  and  not  that  appropriate  to  the  parent  sea  water,  tf.  Further,  Tabata  (Ref.  2)  has  shown  that  the 
chlorinity  and  freezing  temperature  of  the  brine  run-off  is  a  function  of  the  freezing  fraction,  n(=  R,/Rw  ). 
In  fact  his  data  suggests  the  simple  relationship; 

ts  =  (1  +  n)  tf  (B.15) 

although  this  likely  does  not  hold  as  n  -*•  1 . 

The  freezing  temperature  is  related  to  the  salinity  by : 

tf  =  -0.002  -  0.0524S  -  6.00  X  10'5S2  (B.16) 

where  S  is  in  parts  per  thousand. 

This  expression  is  within  0.002°C  of  the  values  given  in  Reference  39  for  salinities  between  0%o  and  40%o . 

Thus,  incorporating  Equations  (B.2)  and  (B.3a)  in  Equation  (B.6): 

qw  =  1.7  X  10’7  Hwcw  V  (td  -  ts)  (B.6b) 

where  td  and  ts  are  given  by  Equations  (B.14),  (B.15)  and  (B.16). 

B.1.3  The  Convective  Heat  Term,  q,. 

The  convective  heat  transfer  is  expressed  in  terms  of  a  heat  transfer  coefficient,  h,  and  the 
temperature  difference  between  the  surface  and  the  surrounding  air: 

qc  =  h  (t,  -  ts)  (B.17) 


In  determining  an  expression  for  the  convective  heat  transfer  coefficient,  a  certain  amount  of 
approximation  is  necessary  since  the  shape  and  size  of  the  icing  surface  is  undefined. 
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For  a  flat  plate  in  turbulent  flow  parallel  to  its  surface  (c.f.  the  deck  of  a  vessel)  an  average 
heat  transfer  coefficient  over  a  length  L  is  given  (Ref.  53)  by: 

ka 

h  =  0.037  —  Prl/3Re08 
L 

If  a  mean  “film”  or  boundary  layer  air  temperature  of 

reduces  to: 

V0.8 

h  =  6.44  -  W/m2K  (B.19) 

L0.2 

In  the  case  of  a  cylinder  of  diameter  D  in  crossflow  at  a  Reynolds  number  of  between  50,000 
and  200,000  the  heat  transfer  coefficient  may  be  expressed  (Ref.  54)  by: 

ka 

h  =  0.0208  —  Re0814  (B.20) 

D 

An  expression  that  provides  the  same  index  to  the  Reynolds  No.  as  in  Equation  (B.18)  and  is 
within  1%  of  Equation  (B.20)  over  its  applicable  Reynolds  number  range  is: 

ka 

h  =  0.0244  —  Re0  8  (B.21) 

which  at  a  mean  film  temperature  of  -  5°C  reduces  to: 

V0.8 

h  =  4.75  -  W/m2  K  (B.22) 

d°.2 

Expressions  (B.20)  through  (B.22)  are  appropriate  to  smooth  cylinders  in  air  flows  of  low 
turbulence.  The  effect  of  roughness  and  of  free-stream  turbulence  can  be  considerable  (Ref.  54),  and 
since  the  air  flow  over  a  ship  is  highly  turbulent,  it  does  not  seem  unreasonable  to  increase  the  coeffi¬ 
cient  of  Equation  (B.22)  arbitrarily  and  equate  it  to  that  of  Equation  (B.19).  If  then  a  characteristic 
dimension  is  arbitrarily  chosen  as  3  m  in  each  case,  there  results  the  expression  for  both  flat  and  cylin¬ 
drical  surfaces: 

h  =  5.17  V0'8  V/m2K  (B.23) 

Because  of  the  power  (0.2)  to  which  the  characteristic  dimension  is  raised,  the  heat  transfer 
coefficient,  h,  is  not  very  sensitive  to  this  dimension,  being  within  10%  of  its  nominal  value  given  by 
Equation  (B.23)  for  values  of  L  from  1.8  m  to  5  m. 

Equation  (B.17)  can  finally  be  written. 

qc  =  5.17  V08  (ta  ts)  (B.24) 


(B.18) 

-  5°C  is  assumed,  this  expression 


where  ta  is  the  air  temperature,  and 

t,  is  the  surface  temperature  defined  by  Equation  ( B.  1 5 ). 
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B.1.4  The  Evaporative  Heat  Term,  qe 


The  evaporative  heat  transfer,  qe ,  can  be  evaluated  by  analogy  with  the  convective  heat 
transfer  (Ref.  52): 


9e 


e  £v 
pcp 


<ea  -  es) 


(B.25) 


where  Pr  = 
Sc  = 
e  = 


P  = 

Cv  = 


Prandtl  number  =  0.711 
Schmidt  number  =  0.595 

ratio  of  molecular  weights  of  water  vapour  and  dry  air  (0.622) 
atmospheric  air  pressure 

latent  heat  of  vaporization  of  water  (2.5  X  106  J/kg) 
specific  heat  of  dry  air  (1.005  X  103  J/kgK) 

saturation  vapour  pressure  of  moist  air  at  temperature  ta,  ts  respectively. 


Evaluation  of 


,  e,  fiv  and  cp  results  in: 


qe  =  1731  ^  (ea  -  es)  (B.26) 

and  if,  as  before,  h  =  5.17  V°  8W/m2K 
and  assuming  p  =  100  kPa 
then: 

qe  =  89.5  (ea  -  es)  V°  8  W/m2  (B.27) 

where  ea  and  es  are  in  kPa,  and  V  is  in  m/s. 

The  saturation  vapour  pressure,  e,  may  be  expressed  in  kPa  as  a  function  of  temperature,  t(°C), 
by: 

e  =  ^((1.9226  X  10'7t  +  2.4545  X  10'5  )t  +  1.4224  X  10'3)t  +  0.044436)t  +  0.61094  (B.28) 

This  4th  order  polynomial  regression  equation  is  within  V*%  of  the  values  given  in  Reference 
39  in  the  temperature  range  -25°C  to  +10°C. 


B.1.5  The  Conductive  Heat  Term,  q. 


In  this  analysis  the  heat  conducted,  qa,  to  or  from  the  icing  surface  through  the  ice  layer  and 
the  underlying  structure  will  be  assumed  negligible,  although  frequently  instances  msy  be  observed  on 
iced  ships  of  inadequately  insulated  crew  spaces,  as  indicated  by  a  lesser  quantity  of  ice  in  these  areas, 
thus  clearly  indicating  that  in  these  instances  heat  is  being  conducted  through  the  structure. 
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B.2.9  The  Total  Heat  Balance  and  Rate  of  Icing 

The  overall  heat  balance  (Eq.  (B.l))  may  now  be  written  using  Equations  (B.5),  (B.6b), 

(B.24)  and  (B.27): 

efRi  +  0.68  HwV(td  -ts)  +  5.17V0  8  |(ta  -  ts)  +  17.3  (ea  -  es)}  =  0  (B.29) 

and  evaluating  C,  as  3.33  X  10s  J/kg,  the  rate  of  icing  is  given  by: 

Ri  =  2.04  X  10-6  HwV(ts  -  td)  +  1.55  X  lO'5  V0-8  j(ts- ta)  +  17.3  (es  -  ea)j  kg/m2 s  (B.30) 
Or,  alternatively  in  terms  of  rate  of  growth  of  ice  thickness,  Nj.- 

Ri 

N,  =  —  =  2.3  X  lO'9  Hw  V(ts  -  td)  +  1.74  X  10'8  V0  8  j(ts  -  ta)  +  17.3  (es  -  ea)|m/s  (B.31a) 
where  the  density  of  ice,  pt,  is  assumed  to  be  890  kg/m3 . 

If  the  rate  is  to  be  expressed  in  the  more  convenient  units  of  mm/h,  then: 

-  8.26  X  10'3  HwV(ts  -  td)  +  6.28  X  10'2  V0  8  { (ts  -  ta )  +  17.3  (es  -  ea)Jmm/h  (B.31b) 

In  these  equations  the  wave  height,  Hw ,  is  in  metres,  the  relative  wind  velocity,  V,  is  in  m/s, 
the  temperatures  are  in  °C  and  the  vapour  pressures  are  in  kPa. 

If,  as  is  more  usual,  the  velocity  is  given  in  knots.  Equation  (B.31b)  becomes: 

Nt*  =  4.25  X  lO'3  Hw  Vkt(ts  -  td)  +  3.69  X  10*2  Vkt08  |(ts  -  ta)  +  17.3  (es  -  ea))mm/h  (B.31c) 

To  test  this  model  of  icing  under  conditions  of  sea  spray,  the  icing  rate  as  calculated  by  Equa¬ 
tion  (B.31c)  using  the  conditions  applicable  to  39  reported  icing  incidents  (Table  VIII)  has  been  com¬ 
pared  in  Figure  43  with  the  mean  icing  rates  derived  from  these  reports. 


B.3.0  Practical  Implementation  of  the  Rate  of  Icing  Equation 

In  evaluating  those  expressions  that  depend  on  the  vapour  pressure  of  the  atmosphere  (i.e. 
Eq.  (B.12),  the  evaporation  factor  in  the  drop  cooling  equation,  and  Eq.  (B.27),  the  evaporative  heat 
term  in  the  heat  balance  at  the  icing  surface),  it  has  been  assumed  that  the  air  is  90%  saturated  rather 
than  saturated,  so  that  in  both  Equations  (B.12)  and  (B.27)  the  term  ed  is  replaced  by  the  expression 
0.9  ea. 


Because,  from  Equation  (B.15),  the  effective  surface  temperature,  ts,  is  a  function  of  the 
freezing  fraction,  n(=R,/Rw  )  and  hence  a  function  of  the,  as  yet  unknown,  icing  rate,  it  is  necessary 
to  employ  a  simple  iterative  procedure  to  determine  the  appropriate  value  of  ts  to  use  in  Equations  (B.6), 
(B.24)  and  (B.28).  The  initial  value  of  the  freezing  fraction,  n,  was  assumed  to  be  zero,  i.e.  ts  =  t,  .  Sub¬ 
sequent  iterations  used  the  computed  value  of  Rs  (Eq.  (B.30))  to  determine  a  new  value  of  n,  until  suc¬ 
cessive  values  of  n  differed  by  less  than  an  arbitrary  amount  (0.0001  was  used).  If  n  assumed  a  negative 
value,  it  was  set  at  zero  and  no  icing  was  inferred . 
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Similarly,  in  determining  the  final  drop  temperature,  td ,  it  was  initially  set  at  the  sea-water 
temperature,  tw ,  and  a  value  of  Xt  (Eq.  (B.12))  computed;  this  was  in  turm  employed  in  Equation  (B.14) 
to  derive  a  new  value  of  td  to  be  used  in  Equation  (B.12).  This  procedure  was  repeated  until  successive 
values  of  td  differed  by  less  than  an  arbitrary  amount  (again  0.0001°C  was  used). 


! 
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FIG.  B.1:  COOLING  OF  2  MM  DIAMETER  WATER  DROPS  FALLING  AT  TERMINAL 
VELOCITY  THROUGH  AIR  OF  TEMPERATURE  -  10°C 
INITIAL  DROP  TEMPERATURE  +2°C 
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